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Zur Menbrnn der Drehschalen negativer :: Krümmung. 


Herrn Prof. Dr..L. Prandtl zum 70. Geburtstag. 


Von W. ‚Flügge in Kresbronn. 


Für Den negativer Krümmung hat die Differentialgleichung der M nenne reelle 
Charakteristiken. Daraus folgen Besonderheiten des Spannungszustandes, durch die sich diese Schalen von. 
den kuppelartigen Schalen positiver Krümmung wesentlich unterscheiden. Sie werden mit mathematisch _ 

h ‚einfachen Mitteln dargestellt, und es wird auf die Notwendigkeit einer genaueren Untersuchung hingewiesen. 
With regard to rotary shells of negative curvature the method of the differential equation of ihe mem- 
. brane intersecting forces has to show real characteristics. hereupon are resulting particularities on the 

 elasticity state by which these shells are essentially differentiating of dome-shaped shells of positive curvature. 
They are produced by simple masbemaztal means; great stress, however, is taid. om the necessity of more 
‚careful researches. \ | 


‚En ce-qui concerne les nike rotatives de cowrbure negative la 'methode de l’öquation differentielle 


des forces d’intersection de membrane dispose de reelles caracteristiques. De la resultent des partieularites et 


de l’etat d’elastieite par lesquelles, ces coquilles en question se distinguent essentiellement des coquilles 
de courbure positive. Elles sont reproduites par de simples moyens math&matiques; toutefois grande valeur 
“est attachte ü la necessitE d’une investigation plus exacte. 


a mepemeHHLX (BPammamınuxca) BEIANBIIeH HETATHBHOTO CRPHB.TEHHA, audbepenH- 

' IHaNBHbIe YPaBHEeHHA COYETATENIBHEIX CHI MEMÖPAHE PACHONATAIMT PeasIbHBIMH XAaPaRkTepuCTH- 
kaMu. Ms 3TOTO BETeKAWT OCODEHHOCTH COCTOAHHA HANPSIKEHHA, SHAUHTEILHO OTIIHYAIOIIUE 

 DTH BEIANBIIIH OT: KYNO000Pa8HEIX" BRITANEINIeH IIOSHTUBHOTO UCKPuBAeHHA. OHn W300pa- 
3KAIOTCH MATEMATHTECKH-IIPOCTBIMH CPENCTBAMN. OEBamAeI BHUMAHHe HA HEODXONHMOCTb 

TOYHOTO HCCJIENOBAHHA. 

In der bautechnischen Literatur ist gelegentlich über den Bau von Kühltürmen aus Eisen- 
beton berichtet worden, die in Form dünner Hyperboloidschalen ausgeführt sind und an ihrem 
unteren Rand auf einer großen Zahl von Einzelstützen ruhen. Aus diesen Berichten ist nicht zu 

erkennen, daß man sich bei der Berechnung dieser Bauwerke der Besonderheiten bewußt geworden. 
ist, durch die sich die Statik solcher Drehschalen negativer Krümmung sehr wesentlich von der 
der Kuppeln unterscheidet. Es ist das Ziel dieser Ausführungen, diese Unterschiede zu zeigen 
und auf die ganze Problematik der Membrantheorie negativ gekrümmter Schalen hinzuweisen. 
Die Membrantheorie der Drehschalen geht bekanntlich aus von den 
drei Gleichgewichtsbedingungen der Schnittkräfte am Schalenelement. 
Diese Schnittkräfte sind die in Schnitten längs der Koordinatenlinien je 
Einheit der Schnittlänge übertragenen inneren Kräfte, nämlich die . 
Meridiänlängskrait N,‚die Ringkraft N, und die Schubkräfte N.5=N 59. 
(Bild 1). Das Gleichgewicht dieser Kräfte und der auf dem Schalenele- 
ment liegenden Lasten in Richtung der Koordinatenlinien und der 
“ Schalennormalen führt auf die drei Gleichgewichtsbedingungen 
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Bild. 


Es sind zwei Differentialgleichungen und eine einfache lineare Glei- 
chung. Mit dieser letzten kann man eine der Längskräfte, etwa die Längskraft N,, eliminieren 
und behält dann ein lineares Differentialgleichungssystem 
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. Randwerte nicht in das Innere der Schale fortpflanzen, sondern sich ausgleichen, \ [ 
von Kuppeln auf Punktstützen bekannt ist!). Hat dagegen die Schale negative Krümmung, 
ist der Koeffizient von 92 N ,/0 92 negativ, die Gleichung also von hyperbolischem Typ. Solche 
Gleichungen besitzen reelle Charakteristiken, d.h. Kurven auf der Schalenfläche, längs derensich 
' Unstetigkeiten der Randwerte über die ganze Schale fortpflanzen. Dadurch treten im Kräftespiel 
der negativ gekrümmten Schale Erscheinungen auf, die der an Schalen positiver Krümmung ge 
schulten Anschauung zunächst überraschend erscheinen und, schon beim Entwurf von Schalen- 
‚ konstruktionen RUROLTISENE erfordern. Aus dem Vorhandensein reeller Charakteristiken ent- 
springt aber auch ein Lösungsverfahren für die Dif- 
ferentialgleichungen. Diese Dinge sollen hierandm 
- einfachsten Beispiel, gem Rotationshyperboloid, vor- ER 
geführt werden. — et 
Die in Bild 2 in rind: und Aufriß dargestellte 
Schale hat die- Gleichung i a 
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Punkt des Kreises z=0 berührt, etwamit derEbene. : 
y=a, so erhält man durch Einsetzen in die Hyper- Ne 
beleiapleichung! für die Schnittkurve er 
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tg ‘= bt „gegen die Grundrißebene. Auf dem Hyperboloid gibt es also zwei Scharen von "geraden 


Linien, seine Erzeugenden. Wir stellen zunächst noch einiges über die Geometrie dieser Linien fest: 
Bild 3 zeigt die beiden Erzeugenden, die durch einen beliebigen Punkt A des Hyperboloids 

gehen. Wir bestimmen den Winkel ß, unter dem sie sich schneiden. Auf der Meridiantangente ist 

eine Strecke AC der Länge 1 abgegrenzt. Ihre Projektion auf die Lotrechte durch A ‚hat die 


2 Länge AB=sing. Aus dem Dreieck ABD folgt AD= un und damit aus dem-Dreieck- RR N 
ADO: Ba | | = 
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Ferner interessiert der Winkel w zwischen den beiden Meridianebenen, die eine Erzeugende im. 
we Kehlkreis und in einem beliebigen Punkt der Schale schneiden. Nach Bild 2 ist 


Für die Krümmungsradien liefert die analytische Geometrie die Formeln 
, a2 b2 a2 
ar (a? sin? 9 — b2 cos? 9)3/2’ 


Yarsin?p —b2 cos? 
!) vgl. z.B. W. Flügge, Statik und Dynamik der Schalen. Berlin 1934, 8.48. 
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Damit kann man auch schreiben: 


gendenschar begrenzter Streifen eingezeichnet. Leitet manin ihn P 
an seinem einen Ende eine Kraft P ein, so ist er und damit die x 
‚ganze ‚Schale offenbar im Gleichgewicht, wenn dieselbe Kraft P 
in entgegengesetzter Richtung an seinem anderen Ende angreift 


“ herangehen. In Bild 4a ist auf der Schale ein durch zwei benachbarte Gerade der gleichen Erzeu- 


und der Streifen dadurch wie ein einfacher Zugstab beansprucht X % % 
wird. Wir haben damit einen Spannungszustand vor uns, der dat Mn 
‚eine Lösung der Differentialgleichungen (1) sein muß, und wir 

können uns auch überzeugen, daß es der einzig mögliche ist, wenn ! AN 
die Kraft P an dem einen Rand als Last gegeben ist und der 4, My, 


andere Rand gestützt ist. Um das einzusehen, wollen wir aus 
der Schale ein Element herausschneiden, das durch je zwei be- 
nachbarte Erzeugende jeder Schar begrenzt ist (Bild 4b). Die in 
seinen Kanten wirkenden inneren Kräfte zerlegen wir in Richtung RES ndehe 
der Erzeugenden in schiefe Längskräfte N;, N, und diezugehörigen 


Schubkräfte N:„—= N,.. Wenn wir. für diese Kräfte das Gleichgewicht in Richtung der. 


Schalennormalen aufstellen, so tragen die schiefen Längskräfte dazu nichts bei, denn sie 
haben ja wegen der Geradlinigkeit der Erzeugenden an'gegenüberliegenden Seiten genau die 
gleiche Wirkungslinie, können also höchstens eine Resultierende in dieser Richtung haben, aber 
nicht normal dazu. Dagegen liefern die Schübkräfte Beiträge, denn die Kräfte N: „ auf zwei gegen- 
überliegenden Seiten sind nicht parallel, da.sie ja die Richtung von zwei V erschiedenen 
Erzeugenden derselben Schar haben, die einen kleinen ‚Winkel miteinander einschließen. Wir 
wollen es uns hier ersparen, diesen Winkel festzustellen und die Gleichung explizit anzuschreiben, 
sondern begnügen uns mit der Feststellung, daß sie die Form haben muß: N;,„ multipliziert mit 
einem aus der Form der Schale folgenden nicht verschwindenden Koeffizienten ist gleich der Normal- 
last Z. Daraus folgt, daß die nur am Rande belastete Schale im Schnittsystem der Erzeugenden 
keine Schubkräfte N;, hat, daß also die Last P, die an einer Seite in den Streifen zwischen zwei 
Erzeugenden eingeleitet wird, nicht seitlich aus ihm herauskann, sondern bis zum anderen Schalen- 
rand durchlaufen muß, wo sie als Einzelkraft auf die Stützung wirkt. \ 
Ein solcher Spannungszustand ist eine physikalische Unmöglichkeit; denn wenn nur längs 
eines schmalen Streifens in der Schale starke Spannungen und daher auch große Dehnungen vor- 
handen wären, so würde er ja gar nicht mehr in den Zusammenhang des Ganzen hineinpassen. 
Da die Membrantheorie nichts anderes liefert, so muß also in diesem Falle die Biegesteifigkeit 
der Schale eine wesentliche Rolle spielen. Wie sich dadurch der Spannungszustand umgestaltet, 
- ob der kraftführende Streifen wenigstens in verschwommener Form durch die ganze Schale hin- 
durch noch erkennbar bleibt oder ob ähnlich wie bei den Kuppeln in einer Randzone schon ein 


weitgehender Ausgleich der Spannungen zustande kommt, das läßt sich ohne Aufstellung der 


vollständigen Biegetheorie nicht vorhersagen. 

Bei anderen negativ gekrümmten Schalen, deren Charakteristiken keine geraden Linien sind, 
sieht der Spannungszustand etwas komplizierter aus. Eine am Rande wirkende Einzelkraft er- 
zeugt dann Schnittkräfte nicht nur längs der beiden durch ihren Angriffspunkt gehenden Charak- 
teristiken, sondern in dem ganzen dazwischen liegenden Gebiet. Die Berechnung wird dadurch 
etwas mühsamer und die Unverträglichkeit der Formänderungen weniger grob. Auf die Einzel- 
heiten der Rechnung, die in der mathematischen Literatur zu finden sind 2), brauchen wir hier 
nicht einzugehen. 

Daß die Membranlösung für eine Einzelkraft am Rande zu Schwierigkeiten im Formände- 
rungszustand der Schale führt, ist kein Hindernis, aus ihr.durch Integration die Lösung für stetige 
Randbelastungen aufzubauen, sobald man dabei auf Spannungszustände kommt, die stetig sind 
und somit die Verträglichkeitsforderung der Formänderung nicht mehr verletzen. Wir zeigen das 
für die Belastung des Kehlkreises mit einer harmonischen lotrechten Last Np = 8, cosnd. Dann 
wirkt auf ein Linienelement a - dd, des Randkreises die lotrechte Kraft 8, acosn® dd, die wir in 
84 I 22.38 in Richtun 

2 sin & 
zerlegen (Bild 5a). 


zwei Kräfte 


g der beiden durch diesen Punkt gehenden Erzeugenden 


2) vgl. z.B. Riemann- Weber: Die Differential- und Integralgleichungen der Mechanik und Physik, 8. Aufl. 


Braunschweig 1930, S. 793 ff, 
5* 
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Nach diesen geometrischen Vorbereitungen können wir an eine einfache statische Überlegung : 
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unkt/aoi n | ie beiden nden, die den 
Durch den Schalenpunkt mit den Koordinaten 9, ® gehen die beiden Erzeugenden, die ( 
kreis ir Meridianen 9 + treffen. Betrachten wir in diesem Schalenpunkt ein dreieckiges Zi 
en 0 n Veelemdut(Bilänp), nase oh Fee n 
gezeichneten ist,so wirkt auf dessen rechte Seite die Bi 
Kraft 
up,z a o)d® 
sin & 
und auf die linke j | 
Br 8,98 n.(0 — o) dd 
2sina 
Auf die waagerechte Seite r,d®. wirkt die Resul- 
tierende dieser beiden Kräfte mit der Meridiankom- 


ponente 
8, add 
N, 0= Sn [eosn (+ ®) 
| ER N 


2 


Bild5a und b. 
und der Schubkomponente 

| .. [cosn (+ w) —cosn (d—o)]sin &. 
Daraus berechnen sich die Schnittkräfte durch Einsetzen der geometrischen Beziehungen Gl. (2), 


(3), 44), ©): 


Ng8 1,0% = 


8a COSNÜCOSnw sin“ COS WCOSN 
N,= N ee EN 
u sin & sin @ sin o 
r, COS WCosn w: a? sin? po —b? cos? 
Na 0 nn  ı 
7, sin 9 b 
8.02 ’ 
Sn 5 cos? wsinpcosn wcosn®, 
N %asinndsinnwo. PB Ss, asinn®sinno Ysin® p —sin? & 
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% 19 sin & 2 ro, sin & sin @ 
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Ganz entsprechend findet man für eine Tangentialbelastung N, s= T, sinn® die beiden Kräfte 


T . ’ ” urE 
Be OL RASTE DEN up AR OnULAN @)dd 
208 «& 2c0S a 
und daraus. für das waagerechte Schnittelement r,d® die Kräfte 
Rad‘ - ß 
N, 09 = Dos, Sinn (+ w) — sinn (d — w)] cos 5° 
The ‘ i 
N, 7,40 = sy Ban (+ ©) rn = 


Drückt man wieder alles durch @, 9 und w aus, so erhält man 


Tb cos @.. 
N,=7- —— 2sinnwcosnd, 
a sin 


a. 
N T,,sinpcos’osinnwcosnd, 


N,» = T,„c0o®?wcosnosinn®. 


Diese Formeln geben die allgemeinste Lösung für eine beliebige Randbelastung. Es ist lehrreich, 
sie mit den bekannten Formeln für Schalen positiver Krümmung) zu vergleichen, z. B. mit der 
geschlossenen Lösung für die Kugelschale mit zwei freien Rändern. Dort kann man nur ein Kräfte- 


®) W. Flügge: Statik und Dynamik der Schalen. Berlin 1934, S.43 bis- 44. 


"ngetirn en BR dem Längskraft Bd Schubkraft für fele reihe in a m Verkaais 8, 
zueinander stehen. Dagegen sind die Kräftesysteme, die man an den beiden Rändern vorschreiben 
"kann, völlig unabhängig voneinander. "Hier, bei der Schale mit, negativer Krümmung, liegen die 
"Verhältnisse durchaus anders. Wir können an einem Rande 9 = 94 beide Randschnittkräfte vor- 
ee 1 geben: N,—=N, a= N,nacosnd, Bere Heinz madı und erhalten dann für 8, eh. 
Ar u, die Gleichungen ER U ER VE RER a EN ge f RL 1 
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8 Bo Diese Gleichungen kann man immer 1oser.d dern) es ist ja immer PA, en 0, 4 Deik 3° ‚so daß die, 4; ei 
ER ee Bar 
N rechten Seiten nicht unendlich. werden können, und auch - die Nennerdeterminante cos? m [Oy IR 


Be + sin?n o, kann nicht O'werden. Über die Schnittkräfte an dem zweiten Schalenrand ist damit 
SL "aber schon verfügt. Dabei zeigt sich im allgemeinen, daß zu einer Belastung des einen Randes A 
7% nur mit eg oder nur mit Schubkräften jedesmal eine : Kombination von Beiuen, 2 
ERENTA am anderen Rande gehört. RE At RE 
Versucht man, ähnlich wie bei den positiv Belruninlen ‘Schalen an een Ba eine ee BSL Bi 
Schnittkräfte vorzugeben, so kommt man,auf ein.anderes Gleichungspaar für $,, Ta, z.B. won 
man an n.beiden: Rändern A und B die Längskräfte geben will: 


‘ f /- N 
” = ee 2 
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a ER das Glichungstyeten so auffassen, daß am- Rande A außer der Längskraft Ngn solche - 
Schubkräfte-N,;,„ anzubringen sind, daß durch die gemeinsame Wirkung dieser beiden Last- 

gruppen am Rande B gerade die vorgeschriebenen Längskräfte N,„ entstehen. Das ist natürlich Ber 
nur dann möglich, wenn zur Schubbelastung des einen Randes eine Belastung des anderen gehört, N 
die auch Längskräfte enthält. Da das zwar in ‚der Regel aber nicht immer der Fall ist, gibt es 
Versager. Die Nennerdeterminante 1 


ge 


b ; A b 
EIERN TO EB AN BAT =, sinn (on —@,) ET  :, 


der beiden Gleichungen kann nämlich null werden, wenn 4 —®z zu zin einem 1 rationalen Ver- 
as steht: 
| N 
Mg u 
Dann verschwindet sinn (v4 —®z) für n= m, und alle ganzzahligen Vielfachen davon, also 
für unendlich viele Harmonische. Der Winkel d ändert sich dann beim Durchlaufen einer Charak- 


-teristik vom einen bis zum anderen Rande gerade um ein ganzzahliges Vielfaches von m. also des R 


Abstandes der Maxima und Minima der nach cos rn ® verteilten äußeren Kräfte. D.h. daß die 
beiden Charakteristiken, die in irgendeinem Punkte des Randes B zusammenlaufen, von zwei 
Punkten des Randes A ausgehen, die zu derselben Phase der Verteilung cos n 9 gehören, also 
Kräfte mitbringen, die sich wieder genau zu einer Meridiankraft zusammensetzen, wenn sie von 
Meridiankräften am ünteren Rande herrühren, und zu einer Schubkraft, wenn am unteren Rande 
nur Schubkräfte angreifen, In diesem Falle ist es also unmöglich, $„ und 7, so zu bestimmen, daß 
an beiden Rändern willkürlich vorgeschriebene Meridianlängskräfte auftreten, sondern diese beiden 
müssen in einem durch die Geometrie der Schale bestimmten Verhältnis zueinander stehen. 


Man kann diese Schwierigkeit zwar formal leicht umgehen, indem man die Schale so formt, 
daß a eine Irrationalzahl ist. Praktisch ist damit natürlich nichts gewonnen, denn in Er 


einer solchen Schale wird die Nennerdeterminante zwar für keine Harmonische null, aber doch für 
manche sehr klein, und dann ergeben sich Werte 8, und T'„, die sehr groß gegenüber den gege- 
benen Randwerten der Meridiankräfte-sind und dementsprechend sehr große Randschubkräfte zur 
Folge haben. Dieser Spannungszustand kann wohl auftreten, wenn man wirklich die vorgeschrie- 
benen Meridiankräfte, die auch an einem Rande nullsein können, aufbringt und der Schale an beiden 


n Bingegangen\ am 8.10. 1944. 
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u A erbiegen, en steif ist! und Sch Be ebig 
i ‚recht stehen. Da aber das Verhältnis der Steifigkeiten nie 


vorziehen, kleine zusätzliche Meridianlängskräfte aus der Schale aufzune n 
» Randschubkräfte wesentlich abgebaut werden können. EIER 


= Um dieses Kräftespiel zwischen Schale und Ring weiter zu verfolgen, | & 
die Formänderungstheorie der Membranschale aufstellen, vielleicht sogar die vollständ 


e 2 theorie für nicht drehsymmetrisch verteilte Randbelastungen. Einfach kann diese we 


denn i in ihr müssen sich j ja die sonderbaren Verhältnisse der Membrantheorie irgendwie auswirk 


Hierauf einzugehen, sei einer späteren ausführlichen Untersuchung vorbehalten. Ziel dieses kurzen 

Beitrages soll es sein, die Probleme zu zeigen, die sich aus der negativen Krümmung. der Schale 
"ergeben, und für das lo Be u das ‚einschalige a die Merabranläsnne zu 
or ‚entwickeln. a ’ 
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_ differentialgleichung ebener laminarer en bei hohen“ 
: | Reynoldsschen Zahlen.) | 

- Von W. Tollmien, z.Z. in Farnborough. 

4: 'Fehlerabschätzung von 93 und 91. 


4,1 Definition und Eigenschaften von o). Das asymptotische Näherungsintegral op, ist nach 
| &. (9) aus BP, Pa und 9% wo) aufgebaut. Um den Anschluß an frühere Untersuchungen v zu 
2 en setzen wir ! 
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| Nach 1 Gl. ( (8) und. 3. Gl, , erhalten wir dann die Diterentigeichung Dr Pa: 


a» dpa TR SEELE 
dn? a dn? U,- ae (2) 
Es läßt sich nun mittels Hankelscher Funktionen [vgl. die folgende Formel (&)]eine Tess 
d? P31ı 
dn?, 
Genauer gilt für große nl: 


von Gl. (2) herstellen, die in der unteren m: -Halbebene sich für große |»7| wie D% verhält 10), 
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pa il an (3) 
| ET ET RE 5 
so daß aus Gl. (2) i 
v d!9 U, 2 
RE dn? BERGE e 
und daraus durch Integration \ 
d’o;; u, 


dr? m Dig $ a R Bo SS N ST N 
folgt. Weiter kann man es einrichten, daß für. große.|n]| 
rose 0% | BE 
A Te ENTER RI a 


e 0 
und : h 


IR 
a ann. 


wird. Danach ist zu beachten, daß &o,, nicht klein mit & wird, sobald |n| groß wird. Die Ver- 


bindung 1 Gl. (9) 
= RB + Bd —1)= Pat EPaı Pı 
geht nun in der Tat bei a [a] in 
1+ E89, = gl) 
N Fortsetzung des Aufsatzes aus Z. angew. Math. Mech. 25/27 (1947), ie 


“) W. Tollmien: a.a.0. 8.27 bis29. In der Anmerkung auf S.29 findet sich ein Vorzeichenfehler, der 
oben in Gl.(3) berichtigt ist. 
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- Über die wichtige aan Gas seien hoch‘ einige Angaben gemacht, Aus der Difterentil- | 
3 Ask. De und dem ‚Verhalten für große nl folgt: at lee: 
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See Del [ entsan ut; Here 
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ur wobei das Argument der Han k e 1 schen. Funktionen er (ü a 12 ist. Diese Lösung wurde von” 


ns mir 1929 qualitativ aufgezeichnet und i in ihrem Verhalten bei großen m diskutiert, 1. H. Schlich- :S 
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- N % ting!%) hat die zweite und erste Ableitung von 9,1 berechnet. Nach seiner Tabelle geben Bild 7 % 


2 N und 8 den Verlauf von a % waere, / Et 
e N 
er 4,2 Reste in der Differentialgleichung beim Einsetzen von Y, und o,. Es soll nun der Rest g, Rr4, 

Et; ‘bestimmt werden, bis auf den 9} die Störungsdifferentialgleichung befriedigt: Bi 
ae N en U —c)(p5" — ap) — U” 9 a 9) .. S 
Dabei müssen die DIRSRESHEIICHURGEN für B,, B, und 95, berücksichtigt Na die wir so SE 
schreiben: 5 8 
RER ER) (U 0) ESTER N ee (10), | 


„ en: U, U— 
19) = (W427 ae 


Dee ea LU Eye ae ee f.urlla)) 


=UR=0. ....dU) 


— pei y=0 gleich U‘ wird. Indem zu Gl. (9) 


| U 
Es ist zu bemerken, daß 


‚iaRL, (Bat ta Rep (Bi) — PBıZaı (921) 
r hinzugefügt wird, erhält man 


i g—4ep, Bı+bep, Pı +4ep, Bi tea HB j 
Di: U, U—c U—c ' | 
—iah—27} PB (U NT U, ER is y E51 Pı 
Brei BUN We UN +). B> 
Bei der Umordnung des Klammerausdruckes setzen wir 
Re en u Se A a Eu n)t 


ı) H. Sehlichting: Amplitudenyverteilung und. Energiebilanz der kleinen Störungen bei der Plattenströ- 
mung. Nachr. d. Ges. d. Wiss. zu Göttingen, Math.-phys. Klasse, Fachgruppe I, Bd. 1, S.47 bis 78, 1935. ‘Vgl. beson- 
ders 8. 57ff. und Tabelle 5. 


Y 
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2 


wobei die Entwicklung von u nach y mit Ir beginnt. Der Ausdruck i in der geschweilten Klammer 


\ 2 
- wird damit ar vi; % a | 
Rumea 24 — 14 
ER RB, (ver. — Roi hge: [72 )+ + yu Pı (eve; U, : a ( s ). 
Durch Integration von Z,, (9,1) erhält man: 


ep, —iaR(U,yep,— age 0 EN Ve NEL Ar TER (15). 


Die Integrationskonstante ist aus dem Verhalten bei sehr großem y bestimmt, wo 9,,, 9,, und 
9, nach Gl. (5), (6) und (7) bekannt sind. Die Beziehungen Gl. RR und (12) werden zur Um- 
formung von Gl. (14) benutzt. Man erhält Sa \ 


9,— EL I DV’ er Pr u B; Ar 4E9, B, 


Kr 6 EP, x an A EP, " i EP; Br: ar Bi t 


Führt man n als Differentiationsvariable von 9,, ein, so bekommt man: 
1 ‚dp u . 1( „d” warn % d? ee) 
et NE, 


d vo 
4% nt € 91 Bir TE ®; 
Der Fehler, mit dem 9, "die Störungsdifferentialgleichung erfüllt, läßt sich sofort hinschreiben: 


oe —iaR{(U —c) Hp) Von Fafti- rn... (4.7): 

4,3 Abgekürzte Darstellungen von 95, 9), 2 und wer Für die Abschätzung: benötigen wir 

noch einige abgekürzte Darstellungen von g,, 9), 9, und 9), die mittels.der Ausdrücke für @,, 
und seine Ableitungen bei großen |n| gewonnen werden. Aus Gl. (3), (4) und (5) ergibt sich sofort 


(16). 


d2 
Eon Be re Na Dee (18a), 
d3 
ni alas tal ee De (18b), 
d!oyı 
RW DREHTE AAN Da NE (186), 
WO 8,, s; und s, beschränkte Funktionen sind. Zur Darstellung von an bei beschränktem % 
führen wir die Funktion z ein: 
T=—Jloge für n)s1, 
1 
t= —loge—1o für 1 —, 
Sen es NR (19). 
te .20 für ml > 
Dann kann man nach Gl. (6) 
doz1 
in Eu HNTEN: HAN) a N er (20) 


setzen, wo t, und t, beschränkte Funktionen sind. Schließlich sieht man an Gl.(7), daß eo, 


bei endlichem y beschränkt bleibt, ohne indes mit e klein zu werden, wie schon.im Anschluß an 
Gl. (7) bemerkt wurde. 


Damit ergibt sich für g, die Darstellung: 


i 1 
= 50 mat eılmhetrimheth BR ae IE N en (21) 
mit beschränkten Funktionen h)...h,. Entsprechend erhalten wir 
9, = Pt €E051 Bı=g Ne Berka 0 lead Suhire Fee ee Segel SE TE A BE NE SE (22a), 
x ‚ dp, c | 
SER an PrtreRı BT ge ee (22b) 


pe 1 d?o,, do , ); 1 
m e dr re EU Arte ST MUurTEtR - „20.2 (226) 
mit den beschränkten Funktionen g,:..g.- 


% 
} 


az Tollmien, Asymptotische Integration. der Störungsdifferentialgleichung 73 

4,4 Schema der Fehlerabsehätzung. Das Schema für die Fehlerabschätzung von Y, und oı 
gewinnen wir durch eine abgewandelte Methode der Variation der Konstanten. Wir bezeichnen 
den Störungsdifferentialausdruck mit 


L(9) = p""—iaR{(U—ı) ("—arpy)—U"P} .. 2... (23). 
In gewissem Sinne kennen wir schon zwei exakte Lösungen gYırı und pıy der Störungsdifferential- 
gleichung; denn wir kennen ihre Näherungen 9, und 9, einschließlich der Fehlerabschätzung. 


Aus den vermuteten Näherungen 9, und 9} und den exakten Lösungen grır und Yry bilden wir 
ein „abgewandeltes‘“ Fundamentalsystem von Lösungen.‘ Es gilt | 


Lo)=g, L)=% Lem=0, Iem=0........ RED, 
Für den Fehler pr —9}5 gilt 
Leu—p)=Lm) —-Lp)=—I---:--- Ber, (25). 
Wir setzen den Fehler mit den variierten Konstanten c,..cy an: 
M—-P=-Hmtaptaputapv - er .. une (26), 
Wie üblich fordern wir. 
VE RR (27), 
(pn) ap top tan tan: rc (27 b), 
(pun—p)"=cPp +%% +6 Pr ee PR ELEN EL 


Aus Gl. (27) und (25) folgen unter Berücksichtigung von Gl. (24) die Bestimmungsgleichungen für 
die ersten Ableitungen der variierten Konstanten: 

‘9 +5% +5 Pım+ ec pn—0, 

9 +5 +5 fmt% pn 

LE HE + Pt nit 

aM +5% +% Pu tPn = 97a a9 
Nur die letzte Gleichung unterscheidet sich durch das Auftreten von c, und c, von den Glei- 


chungen in der üblichen Methode der Variation der Konstanten. Wir setzen die Wrons ki- 
sche Determinante von 9;, 9), Yım und Qıv: 


(28). 


Dr Ps, Yıın IV 

| a 

P "P, Yın Pıv 

| m. m‘ Pin iv 

Ferner führen wir Abkürzungen für Unterdeterminanten der Wronskischen Deter- 
minante ein: 


9% Ym Pıv / Pad; Min, SAN 
DEN en la. De AP. Ayla (806), 
% Pin Pıv Ay, Pa Pıy 
Yı 9%  Pwv Pi Pi Pi 
D=)M m Pr... (80e), D,=1P:% . Pr. a2 (80): 
PM  Piv 4 PM Pu 
Dann ist 
; D 
‘= (gtant% 9) SR a PR NE (31a), 
D | 
= + ıNFt %I) er RER at alle (31b), 
: D 
‘er ra h tt 9) DI EN ER BD ee (Ste); 


D | 
= H+aHNF% 9) nn ER OL (31d). 


. 
» 
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Für die Wronskische Determinante leiten wir mit einer klassischen Schlußweise eine Differential- 
gleichung ab. Schreibt man formal { | ‚N 
L(p)= pP" +a,p"+a,p"+asp +asPp, 


so kann man die vier Beziehungen Gl. (24) als lineare Bestimmungsgleichungen für a)... a, auf- 
fassen. Man erhält speziell 
d W . D 
dy D D Ne a 
Mer er re | ; 
Da a, nach der Definition Gl. (23) von L verschwindet, so ist 
dW 
ee Be ee ae Km hr REN Ari RE (32) 


2455 Darstellung der Determinanten D,. Da. D; und D, u Ausreehnung von W. Nach 3 
G1. (39), 3 GI. (44), 3 Gl. (46) und (22) wird 


£ 17 
Di=750-1(A) 


Zn (H nn 


; 
De ae 


A [Pin Pıya — PıvıPına] 
..(33 2). 


En (792 + 95) [Pın Prva — Pıy Prise] 


1 
2 E o,MutrrTgs+ 7 [Pr Pıvı —Pıy Pr] 


Entsprechend wird, wenn man noch 1 Gl.(4) zu Hilfe nimmt, 
ZEN 


#9 
= ei PBı [PımPıv. — PıvıPrns] 


ar 9[lPınPıyv —Pıypur]) % -.. 2.» (34a). 


| 

a 

a 

& 

ER 
: 


—c 
i DE : 
+70-(8) m: 9, [Pr Pıvı — PıvPımı] 


“ Weiter wird 


1 uU’y N1la 2 
D,= exp Er + 2a Br ca en Pıval9ı (7924 95) — 9,91] 


1 UyY \3e 1 hr 
rt (A) | Pıvı Me HM AurTEt+ «) RE ai . . (89a), 


+ 0-51 (4) | r N pe B ( -,(M)4+745+ a) 9, + 2) 


| 2 1 u,r 
D,—exp ae +) a ou = : 


1/4 
= A Piz [91 (792495) —Pı9ıl 


1 D,ENae 1 En ? 
Ing 0_3/4(H) 7) Bun Pi BE MATTE H+ 2) 94 A . (36a). 


+ 05/4 (H) (ei . co, N)u+rTG+ «) 9 (79%+ a) 


Diese Darstellungen wollen wir an Hand der Tatsache vereinfachen, daß wir die Abschätzung in 
dem endlichen y-Intervall y,...2 Y„ vornehmen wollen. Dann ist 


»al<h<alAl, *»>0 „>0, 
und man kann in den erwähnten Darstellungen o,, (n) durch o, (H) ersetzen. Überhöhen wir 7 


a 1 
durch die obere Schranke log 2,80 kommt man leicht zu den einfacheren Darstellungen 


ng r e) 

Cr 

-i 1 2 
Fn 


Fir 


wo in a beschränkte Funktionen ind A ER Rt mE re 
EN Die Wronskische Determinante 12 soll mittels der Ditferentialbezichung Gl. (32) undeines 
Se Anfangswertes W, bei y= 4, berechnet werden. Dort möge 9% schon praktisch gleich 9, sein, IR 
. ebenso mögen die ersten Ableitungen Eu Dies a nach Gl. (3) bis auf einen Fehler 
«sder rap aorünnne & wenn bei y=Yy. i ER 


nen Ole 


ist, was dureh die bereits Bankeie Voraussetzung. 3.Gl. (6) erfüllt wird. Tanker der schwachen 
Bedingung 3 ‚Gl. (6) können die höheren Ableitungen von 9, PA Ya noch — eo werden, ® TE 
wie en, p, wie &2B. Trotzdem kann man 


4 = £ f f 3 nr j 
Be a ya Pi Pıv 
Al SE Ber Y War 


pP 


im: Pl, 
KR 9, 


‚ ‚ 


Pr Pıy 


setzen. "Während sich dieser Näherungswert sogleich als von der Größenordnung EI? herausstellen 
wird, so sind die vernachlässigten Glieder nur von der Größenordnung &, wie man leicht aus 


‘aus den Ansätzen 1 Gl. (4) und‘1 Gl. (5) bei kleinem FI Unter Ba von 3 Gl. (33) an Stelle 
. der zweiten Determinante | in Gl. (37) ‚ergibt siehe" 


der Bemerkung über p%', X und aus den Darstellungen 3 GI. (39), 3 Gl. (44) und 3 Gl. (46). ir 
für Orr, Yıv und deren Ableitungen erkennt. Bei der Berechnung des Näherungswertes Gl. (37) 


gehen wir von der Beobachtung aus, daß die Wronskische Determinante von 9, und 9, nach 
(der Differentialgleichung 1 Gl. (3) für y > 0 konstant ist. Ihr Wert —1 ergibt sich am leichtesten 


4 * z 5 LE? 


Die Verwendung von 3 GI. (33) bringt dabei Fehler mit sich, die höchstens von der Größenordnung 


e-* sind, wie man aus den abgeschätzten Fehlern von' ‘0, 9/, 9), 9, und @/” schließen kann. 


‘Die Änderung von W im Intervall y,... 2, bestimmen wir durch Integration von Gl. ee 
Diese Aeeation wollen wir mit der a H vornehmen, so daß 


'W= nn \v D, Fanze| 


ist. Für 92, D,, D, benutzen wir "ie Darstellungen Gl. ER (33), (34). Die Abänderung von W 
gegenüber W, ist dern nach Gl. ( at leicht als höchstens von der Größenordnung e”? abzuschätzen, 


an 


92D 


so daß man 


setzen kann mit einem Fehler, der nach den Erwägungen bei der Aufstellung der Formel Gl. (38) 


für W, höchstens von der Größenordnung ®% ist. 


4,6 Fehlerabschätzung von p, und 9}. Für die Fehlerabschätzung von 9}, die wir im Inter- 
vall y,...2y, vornehmen wollen, müssen wir c,, ,, c, und c, nach Gl. (31) ermitteln. c, und c, 
nehmen wir aus Gl. (31a) und (31 b), indem wir als Anfangsbedingungen ce, = c;, = 0 bei y= %, 


- vorschreiben. ‘Die erste Näherung für c, und c, gewinnen wir aus den verkürzten Gleichungen: 


nr D 2 9. D 
| vr a I — ai RE ER (41) 
Nach GI. (21), (33), (34) und (40) erhält man, indem man sich die Integration wieder über H aus- 
geführt denkt: 
B €] =E O4; Ca = &2Q, su Phte, \err se, (wo N Eee. (42) 


mit den; beschränkten Funktionen @, und Q,. Wenn man diese erste Näherung der Lösungen 
von Gl. (31a) und (31 b) nach der üblichen Methode der sukzessiven Approximation verbessert, 


* erkennt man, daß an dem Ergebnis Gl. (42) sich größenordnungsmäßig nichts ändert. Wir merken 


noch an, daß Q, im Intervall %,...2%, mit e*® klein wird. Dies läßt sich zeigen, wenn man 


H D, il, Bunte exp Jah a au 9’ = A,(H) - | h y Es a 65), R 


BR 2 EUR PN {er tn , Re | REN 
ee de 1 sul a.am en oe 


ee 


- 
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bei den Anschatziaden des Rn das man aus Gl. (41) für c, bekommt, die On unter der Be- 
dingung 3 Gl. (6) nach ihrem definitionsgemäßen Verlauf berücksichtigt. 
c, wird sodann aus Gl (31c) und (35) mit der Anfangsbedingung c; = 0 bei y= y, mittels 


der Enstmel . Gl. (28) abgeschätzt. Es ergibt. sich 
v2 ne] 22100 8 
63 —eXPp +89 log 73. EDER REDEN (43). 
Entsprechend wird aus Gl. (31 d) und (36) mit der Anfangsbedingung c,—0 bei y=2 y, 
Ca exp ehe (1+ De log r, DE (44). 


r, und r, sind beschränkte Funktionen von H. 
Nach Verbindung von Gl. (43). mit 3 Gl. (39) und 3 Gl. (44), von Gl. (44) mit 3 Gl. (46) wird 
somit aus ge] Fehleransatz Gl. (26): 


h 1 
91 = EQıPıt 9,9, + e?log 3 04 SR FERNEN Kran (45), 


e * ’ * 1 
Mm —p = EeHNtTPP +E log EEE (46). 


Alle Q sind beschränkt. Im Intervall y,...2 y, werden, wie man aus einer ähnlich schon oben 
angedeuteten genaueren Abschätzung unter der Bedingung 3 Gl. (6) erkennt, @, wie e?/?, Q, wie 
e’/3 klein. 


4,7 Fehlerabschätzung von p, und p,. Für den Fehler von 9, machen wir, wie für den von 
op}, einen Ansatz mit variierten Benin. 


9 M= 619, + 94 039111 +&4P9y ’ 


wobei in Analogie zu Gl. (31) sich jetzt für die variierten Konstanten die folgenden Gleichungen 
ergeben: 


b D 
= (9ı+C19ı+ 629s) mw’ 
D 
1% „= —(91+ 6191+ 639;) nz 
DSH RENNER EN (47). 
= (Jı+ 61914 0292) 7» 
D 
— (9146191 0292) = 


Den variierten Konstanten werden dieselben Anfangsbedingungen wie in 4,6 auferlegt. Es ergibt 
sich dann ganz entsprechend: 


i 
c=e?log m Senne (MB a N Er gi (49), 
exp ie (143) zen] OB (50), 
€ 
a=erp| 1, Y2ar+ we] &tlog 28, TE NE (51). 
Nach Verbindung von Gl. (50) mit 3 Gl. (39) und 3. Gl. (44), von Gl.(51) mit 3 Gl. (46) erhält man 
1 L El 
91 —9P1= e3log = 819ı+ 88,9, + etlog = NER NN RER le (52), 
ARE: I 218, w 1 
9 —9,—=erlog Zap + 828,95 + e?log in SED TEE (83); 


Alle S8 sind beschränkt. 


T 
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5. Zusammenstellung der Fehlerformeln für die vier Näherungsintegrale 
und ihre ersten Ableitungen. 


Damit ist das Hauptziel dieser Abhandlung erreicht, und wir stellen die wichtigsten Fehler- 
formeln zusammen. 


Nach 4 GI. (52) und 4 Gl. (53) ist im Intervall y,...2 Ya 
1 2 EEE 
Mr WE fi er etlog — 85 N EST WER DAR (1), 


1 
9, —9, = log — 281g, resR + log 8, EL Nr er 


wo'84,..19% beschrfrkte Funktionen sind. 
Nach 4 Gl. (45) und 4 Gl.\46) ist in demselben Intervall 


Yu —p, = ee 20,9, + log— % SEHR Kr Nun RE (3), 


uhr = Qi + Ep + E18 9. RR RN Dil as, 


Qı..Q, sind beschränkte Eunkronen, ‚Speziell werden im Intervall y,...2 %, 


0, 03) 9, 0 Ole (5), 
falls dort nach 3 Gl. (6) EN 


Y2> us 
ist. Wie vor 4 Gl. (37) bemerkt wurde, ist dann in dem Intervall %,...2 Ya 
DEE (Eee ERS She EEE ur NR SR (6), 
| WE EONERER NE ea Te (7). 
Für 9, und o, benutzen wir nach 3 Gl. (27) und 3 GI. (26) die Darstellungen: 
IRRR 2 A 
9 = oT exp er (1-+) zon| p; A Ze en NN (8). 
1 vr\su a 
Dt (ep ao zer] BER 9). 
Nach 3 GI. (39) und 3 Gl..(44) hat man. ferner die Darstellungen: 
UND y2 > 
YI = O_5/4 (AH) exp [= (1 u) ne] Pııl (H) ARTEN: (10), 
1 Ur 2194 2 
Yin = ONE (A) Do eXp Bir 1+ t) m Pıtlı (H) Ele a (11). 
Im Intervall y„,...2%. gilt dann nach 3 Gl. (43) und 3 Gl. (45a): 
eP; |P3Im 
IPııı —P3 |< 1 Ep a ee Pe Ba RE te a (12), 
2.221 Ps 
[Prırı — Ede |< u he LS LEN DT, Kae (13) 
mit zwei von e unabhängigen Konstanten P, und P;,.. 
Entsprechend hat'man nach 3 Gl. (29), 3 Gl. a und 3 Gl. (46) die Darstellungen: 
u,r.\14 3/2 
0 = ar4lH) U) @P ir HM DE A, ee: (14), 
1 U,Y \31 2 
9 = 0 (H) m exp ik (1-+) | Ba LER ae air ee (15), 


Dry \oa 
gV= 0254 al exp Bar ya| PvE) ne (16), 


3 
y \3l4 2 
° ep | a 1a O2 Se NR (17). 


j a se, 1 es ven y= oo, Bei i Grenzschichtprefile 
| am (0) bei großen "Wandabständen ee also bei y=® 
Grunde kann in der gesuchten linearen Verbindung ‚der vie sympto 
_ vorkommen, worauf wir schon in 3, ‚3 im Anschluß an die Formel 3 Gi 
Erfüllung, der Randbedingung bei Ur =@© er stellen Br de Sr 
Sr im era) Ya - a > Ya auf. Re as e wE 


er 


SE ist U nach 5 Gl. OB 
yE , 378 o ae, Bir ar: % >” ur : ee I he « 
Sr genügt der reibungslosen Storingsdittsenfialiechugg 1 Gl. (3), läßt sich Also in ee bezei 
neten Intervall aus Por! und Poz» 3 dort ein RT. von 1 Gl. (2) bilden, Eee Rn 
IR; ‚setzen: RR t NR RR Eh Bach ae % 2 ERS 
Ki DIN N  meruBa true Teatane ALLE RT Be re ar. 
wo Yıı Erd 7 nella sind. Nach 1,3 muß Yin dem Intervall RSS ya durch . eine lineare 
Verbindung von 9,1, 93 Psg und 94 darstellbar sein. Da diese vier Funktionen ganz verschie- 
dene Wachstumseigenschaften in dem Intervall ‚haben, so muß jedes der vier Korrekturglieder, 
die nach Gl. (1) zu 9, hinzugefügt werden müssen, um p zu ‚bekommen. i in GE ganz en Inter- . 


vall wie &° 1, klein werden Im), Also erhält man die Darstellung - SR. EZ ee e 2 5 . 


Yu Pt Yır + og Rn —(e ut & Put & Fat Ei =) r ie ER RE nn 3 


Dabei sind &11, &ı» &ı3 und &14 Konstante, Se Inindestene mil £ klein “Werden. Der Nenner 9, (u 


en. EB, | ap den Wert von o,, in der Mitte des Tea also bei = Ya bedeuten. Er mußte wegen des 


überaus raschen Wachstums von Psa hinzugefügt werden. RS Ne Br x x EE: 
Nach 5 Gl. (3), 5 Gl. (5) und 5 Gl. (6) ist im Gebiet y....2 Ye 


Ppu=Pp+O(E) RE a £ Sn Ehe ra SER: (A), a En: 
Genau ı so wie oben könden wir mit, zwei Konstanten y,, und ya SE TEE 


Be Bit : P2 = Yan Par + Yon Pas EEE EEE EL EEN K (5) A 
A setzen. Schließlich ergibt sich BL 


Ar Yzı Ps1 si Ya2Ps2 + &ı Psı a ei, Ps2 + ee, Ps3 & Eu Fr a IS Fa Se Are (6). er EN 
se Le t } 
| wo die Konstanten ar E29, 95 Und &;, höchstens die Größenordnung e haben.  - e 


Nach 5.G1. (8), 5 Gl. (10) und 5 GI, (12) ist im Intervall Ya...2 Ya: ER 


U,Y\%4 EEE BEE 
An Pat 0-,.(H () ep" +4) Hs] 016) BERNER EIES 
Dafür kann man nach 3 Gl. (8) schreiben: | | 


\ Pr Pat d3 9,30 (e) 2 BR de ie N 
Nach 3 @l. (9) ist Bar yo 


= d9;llt oe) er eo 


R din dspalit ot] on. nl 


'*) Umtändlicher kann man diesen Schluß begründen, indem man das 


Also ist 


Korrekturglied of: log .)r von (1) durch . 


eine lineare Verbindung von Ps Ps. 95, und 9, in vier Punkten des Intervalls Ya . 


.2 Ya ersetzt denkt. Dann hat 
man das lineare Gleichungssystem für die vier Konstanten der lineare \ 


n Verbindung zu diskutieren. ° 


_ 


Re gr He „= Se + ’® er Pa Zi Pr + a m Fu | h ur = Bye 
E = a 34 RE a Meei-jet 


es ei at: 
er fi Her Wert von gut Se > U Em Ep En und. Ess sind Konstante, De mindestens wie € klein Sa x 
RE reden, u258 = 

% . Nach a al. 0) ist i im ganzen untersuchten Intervall Ya. v .2; Ya : 


NE 


As DR ER 
Se Det) ("e exp & A Ir 
BEE mit besehränktem Pıy- Daraus ergibt sich für Ye =Y = Ya: Re 
= nn ee 51a u In . 1 R u- d ee 
ee N Ave Ver] Fr [net ie aan ke rt y nenn ae (12). \ Ya 
e < ur %: . k N) . M 


ES e 5 Datür kann ı ‚man ı nach 3°C), a schreiben: 


. 3 R . - = h = & r a“ ya e24 r "up Pıy Er ER . BE £ 2 E 5 . . EN 
5 E ge will man eine Darstellung von Yıv als Tape Verbindung von Ds: 1» 952 Ps3 und P;4, SO Hatı man 
schon in e, nach 3 Gl. (50) und 3 Gl. (51) den Faktor von 9,, bis auf einen relativen Fehler der 
Größenordnung €. Da pıy beschränkt bleibt, so können nach Gl. (13) die Bestandteile mit 0; 


Be 9,5 und 9,; nur die Größenordnung e& Pas „ haben. Also kommt man im Intervall Yar».2Y. ZU Er FE 
ET ‚der Darstellung: = 


A; ER ee 


wo die RES Ey mit e klein wird, während die Konstanten Q,, a, und a, bei beliebig . abneh. 
mendem & beschränkt bleiben. x 


6,3 Korrektur von Pr zugunsten der Randhedingung bei y= x. Der Randbedingung bei En. 
y= & wird 9111 beinahe gerecht, da nach Gl. (11) sein hauptsächlicher Bestandteil das abklin- 
gende 9,, ist. Die hinderlichen nicht abklingenden Bestandteile zb 9,5 und 94 Sn? wir nun 
beseitigen, indem wir eine Korrektur 


Ay dp + OP EEE 1] 


von pm abziehen. Dies führt nach Gl. (11), (6) und (14) zu den ‚beiden Gleichungen für die Kon- 
- stanten Ög vet Öga: 


d; a Ög2 (Yaat 82.) + Öz4 644g Dga ee (16a), 
N € % : s 
a ob ee ee . (166). 
sa ' 
Man erhält 
; 6,2 IR ua NE (17a), 
V22 
d BoszE 
(0) ul ae = 2 Re Re EEN 17b 
= 64 Ps E Y22 


Falls Yos etwa 0 sein sollte,.was in den Beispielen allerdings niemals der Fall ist, so hätte man 
Ö5, p, an Stelle von ö,, 9,, zum Aufbau von A 9,7), zu wählen. 


Yı —A Yin wird im Intervall y,...2 %, gemäß seiner Erzeugung proportional zu einer 
‚ linearen Verbindung von @,; und 9,,, so daß diese exakte Lösung $y1r —4 Yıır über y„ hinaus 
durch eine im Unendlichen abklingende Lösung fortgesetzt wird. Im Intervall y,...%, wird 
dagegen A 977, mit d, 9,, e klein, wie aus Gl. (17) folgt. d, wird unter der Bedingung 3 Gl. (6) wie 


_ (a R)-N2 exp ae z klein. A p/,, wird im Intervall y,... Ya wie d; 9,58 1og e klein. Dem- 


- nach ist es im Intervall %,,...%, für die Fehlerabschätzung von 9, und @, belanglos, ob man 
von der exakten Lösung oder von der korrigierten exakten Lösung 9ırı 4 Yııı ausgeht, 
die im Unendlichen abklingt. 


hi 
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6,4 Zusammensetzung und Korrektur von Yı und gr, um Abklingen bei y= © zu erzielen. 
Aus $ı und Qır soll eine. lineare Verbindung Er \ 3 
| | ort Domina nee EN EN" ER 


gebildet werden, die im wesentlichen einer abklingenden Lösung entspricht. Die nicht abklingen- 
den Bestandteile sollen dann durch eine Korrektur an gr und 9ır beseitigt werden. Rs 


Denkt man sich pı im Intervall y,...2 ys durch @,, Yır durch @, angenähert, so wird 


0,914 0,91 > C,9ı+ %p= OıYıPfatYıaP2)+ C;(YgıPsıt Ya2Ps2) - - - (19). 
Setzt man abkürzend bei y = %,: 


PaataPpa= Da - een a: 9 
at pad nn. ee 2200) 

und nimmt dort als Randbedingung Se 
OD HD IE A Tee (21), 
so ist dies gleichbedeutend mit er 
Ola Bay en ee ee (22). 


\ . 
' Danach bleiben in der Verbindung Gl. (18) an Bestandteilen, die nicht für „=: oo nach 0 gehen, 
noch übrig: 


1 3 Eu 2, Eu 
C,log — (dep: 2Y En 0, ( + en Eee (23), 
E Pas Ps 
wie man aus Gl. (3) und (6) erkennt, Von or soll 
A 9— Ö15911 + Ö14 91V are, Re ker Kim im nl Ina) rel Parka te (24) 
und von 97, Soll 
A 91 Ög9rı + Öza Diva Kerala) Jaltzen atmen Ip ee ee (25) 


abgezogen werden, so daß die nicht abklingenden Glieder Gl. (23) ausgemerzt werden. Diese 
Forderung führt nach Gl. (6) und (14) zu. den Gleichungen: | 


1 5 NER 
log AR Ss dralYıat &)t b10%99 Pa > ren (26a), 
WEN (26b) 
z Des 2, 1464 Ball Me ee ee ee > 
Oo lYaat Ei) Os ala. a ren (27a), 
2% Eu 
== = 69 + Ö6g5a(1-+ &ua) at a Te een (27b). 
Psa Ps4 
Aus Gl. (26) erhält man 
log } 
€ 3 3 
be — (U, —,H) ee. (28a), 
22 
log = 
Öl BF I—— (Ya HT EA): nee. 
u, Yon Por Ya2 &) (28b). 


Im Intervall 4,...y. unterscheidet sich danach @, — Ay, von or bis auf einen beschränkten 
Faktor um e?lög ii 95 9 — Ag; von ; entsprechend um e?log-—p‘’. Daher ändert sich an 
€ 


der Fehlerabschätzung von 9, und o\ nichts Wesentliches, ob man nun gr oder 97—Agral k 
Lösung zuordnet.. Aus. Gl. (27) erhält man A ji ES 


2 2 


da ee N (29), 
2 2 
ds ET A 
ey Fa (29 b) 


Im Intervall %,... y„ unterscheidet sich demnach 9,7 —A Yr von Q7, bis auf einen beschränkten 


Faktor um e? 0% entsprechend 97 —4 Yıı von pr um &2 9%’, so daß sich durch die Korrektur von 
Pır wieder nichts Wesentliches an der Fehlerabschätzung von 9} und % ändert. 


Kl 25 
> 4.071 
Men 


Brit Atıe den. GL: ER folgt 


6, 5 ern an er er N achelein rn ste daß Yyır Be ip Pı gi 0, Yır durch“ 
an die i im Intervall y,...y. nichts ausmachen, zur vollen Befriedigung der Rand- 
 bedingung bei y—= oo gebracht werden Innen, bleiben nur die Randbedingungen an der Wand 
(y= Y,) noch zu erfüllen. Dort muB 9 = —=0 sein, da Han a an der 
Wand verschwinden müssen. Also ist bei y= y, 


x 


ES j = | N 1 Rs Pier ck han n0 RN A Ko 


Out Ort Os Yrn = DR rg (30b) 


Y "zu fordern, wozu die schon aufgestellte Randbedingung Gl. (21) 


ED OD N 


ER ET ’ 1; PET | 
Dr Pin Pr Dre) = Pia Pırw— PraPın EB RN 


_ Diese komplexe Gleichung gibt den Zusammenhang zwischen den drei Parametern R, a, c bzw. 


&, &, c. Sie soll unter Vernachlässigung von Gliedern, die stärker als e verschwinden, angenähert 
Ben, indem 91, Prr und gYırı durch ihre Näherungen ersetzt werden. Im Rahmen des Bee 
Degen wird nach 5 Gl. (1), 5 Gl. (2), 5 Gl. (3) und 5 Gl. (4): ; 


Pıllw,__, Psw un Wosle 
GE En Ne NR ER HR 52), 
Pur w es | “ 
dn/n 
Yu Pin 9 wr zo N a ee: ED FERERRNE 1a (33) 
Y1w”T Bo EeQıw ke nF Dr 2 2% Nr eat ie (34) j 


_ Damit wird aus Gl. (31) 


€ en Bo D.aPıut EeD,.RıwP „) = DraPyu— PzaYın + ID Bar ler 
5 ") 

Die Diskussion von Gl. (35) beschränken wir auf Profile ohne Wendepunkt (U’’ <0). Dann 
stellt sich heraus, daß die Parameter c und & bei R> »bzw. e—>0 so miteinander gekoppelt 
sind, daß ®,„ und 9, nach. O0 gehen. Dies geht erstens aus den zahlenmäßig durchgeführten 
Stabilitätsuntersuchungen hervor. Zweitens wird dies durch einen Satz nahegelegt, den ich bei 
vollständiger Reibungslosigkeit (e = 0) bewiesen habe PB), wonach bei Profilen ohne Wendepunkt 
das gesuchte 9 gleich 9, sein muß, was das Verschwinden von ®,, und 9,,„ fordert. Daher geht 
in Gl. (35) das letzte Glied e®,.„ Qıw Yıw Stärker als e nach 0. Wir können also mit. dem ge- 
wünschten Näherungsgrad die Gleichung für die Bestimmung der Parameter so schreiben: 


€ ER (D..95,—DaP, ED De Da Bin Sr ee (36). 
) 


Näher soll auf das Randwertproblem in dieser Abhandlung nicht eingegangen werden. 


7. Die Störungsdifferentialgleichung bei Grenzschichtprofilen. 


Es war eingangs in 1.1 erwähnt worden, daß U im allgemeinen keineswegs genau eine Funk- 
tion von y allein ist, wie wir doch bei der Störungsdifferentialgleichung annahmen. Schon 1929 
- (a.a.0. S.40) hatte ich indessen andeutungsweise darauf hingewiesen, daß bei Grenzschicht- 
_ profilen der Einfluß von x auf die Grundgeschwindigkeit die Störungsdifferentialgleichung nicht 
wesentlich abändert. Diese älteren Überlegungen sollen jetzt in erneuerter und verbesserter Form 
wiedergegeben werden 1%). 

Die Grenzschichtströmung, deren Stabilität untersucht werden soll, habe die Stromfunktion 
Y (x, y); sei die Bogenlänge der beströmten Wand, und y sei normal zur Wand. Mit 
y (x, y,t) werde die Stromfunktion der Störung bezeichnet. Wenn der Einfluß der Krümmung 
der Koordinatenlinien vernachlässigt werden kann, was wir in Analogie zum Vorgehen in der 


22) W.Tollmien: Ein allgemeines Kriterium der-Instabilität laminarer Geschwindigkeitsverteilungen. Nachr. 
d. Ges. d. Wiss. zu Göttingen, Math.-phys. Klass, Fachgruppe I, Bd.1, S.80—114, 1935. Vgl. besonders $5 und $7. 

1) Zu dieser Darstellung wurde ich durch eine Anfrage von Herrn G. Hamel vom 24.3.1943 angeregt. 
Wesentliche Teile meines Antwortbriefes werden oben benutzt. 
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Grenzschichttheorie einmal ohne Diskussion hinnehmen wollen, so gilt für die N Strom-: 
funktion y= X+ y die Differentialgleichung: 
oAy YyoAy ypoaAy A 
"at "2y Yy dx 900 Yy 
v ist die kinematische Zähigkeit. Glieder, die quadratisch in der Stromfunktion der 
‚ Störung sind, sollen vernachlässigt werden. Dann ergibt sich: 


REDE RI AN 


' ay dx ar 9y | 0 DR 0). 
I A or RA NEN oPoAy 2y ar As 
öy 0x 0y 0% 9c 0% 0x dy 
Die erste Zeile, die sich allein auf die Grundströmung bezieht, verschwindet für sich. Die zweite 
Zeile liefert die Störungsdifferentialgleichung. Bis hierher ist noch alles allgemein gültig für ebene 
Strömungen ohne merklichen Einfluß der Wandkrümmung. 
Jetzt sei daran erinnert, daß bei Grenzschichtströmungen X weit stärker von y als von & 
abhängt. Entsprechend wird die Störungsstromfunktion y hauptsächlich durch die‘ Veränder- 
lichkeit der Grundströmung mit y, weit weniger durch die mit x beeinflußt. Diese qualitativen 
Aussagen sollen nun formelmäßig gefaßt und-in die allgemeine Störungsdifferentialgleichung 
Ay 0oWoaAy yoadAY 2WPoAy 9yoaAyP 


N Bee = LE ORTS 
Bra DRS Pr are et (8) 


einbezogen werden. 

Es werden Dimensionslose eingeführt. Die Geschwindigkeiten werden auf eine charakte- 
ristische Geschwindigkeit U der x-Komponente der Grundströmung bezogen, wobei man etwa 
die Geschwindigkeit nehmen kann, mit der der ruhende Körper vom Unendlichen aus angeströmt 
wird. Soweit .es sich um die Störungsstromfunktion handelt, machen wir alle Längen durch eine ı 
charakteristische Querabmessung der Grundgeschwindigkeitsprofile dimensionslos, die mit der 
Grenzschichtdicke vergleichbar ist. Die nähere Definition dieser Bezugslänge 5 schieben wir auf. 
Auf diese Weise hat man für die Dimensionslosen mit dem Index 1: 


b $ 
wa, d, nn y=lD, I=17 DEUDDR U ee (4). 


Bei der Grundströmung beziehen wir dagegen gemäß der Grenzschichttheorie x auf eine 
charakteristische Abmessung / des umströmten Körpers (z.B. Radius beim Kreiszylinder), wäh- 
rend y wieder auf 5 bezogen wird. Für 5, das eine charakteristische Querabmessung der Grenz- 
schicht sein soll, kann man dann etwa wählen: 


v 

Die anderen möglichen Festsetzungen von 5 unterscheiden sich von dieser nur um einen kon- 

stanten Zahlenfaktor. Es werden also die Dimensionslosen mit dem Index 2 eingeführt durch: 
sen N yayb—ydiN Zaaumbr, (0 Var as (6). 

Dann gilt der grundlegende Ähnlichkeitssatz für Re ai daß 9, (x,, y,) von der Rey- 

noldsschen Zahl unabhängig ist. Für y, (2, Y,,t,) gilt natürlich nichts Entsprechendes. 


Werden diese Dimensionslosen in die allgemeine Störungsgleichung (3) eingeführt, so erhält 
man: 


ae, 9 Ay in (TE boY,9A,yı 
ot, Ye 0%, oy,l® dx, löx,oyl 10%, °Y, 7) 


el b° re 34 
9x, 1? 922 09, v5 191 
mit 
0° 0? 
AR 
tag 
Umgeordnet wird daraus mit: 
oAıyı , Pr 9 Aıyı ap P, > hj Ay oo, 9A,ıYyı 
ot, 9y, 98, 9x, 0y5 Ubl! a a 9%ı 
or OP 0 ps: 9A, ONE OR ln RL > De 8). 
0m,9yy0yı Ubaxloy,ox, U? 9x3 ie 


2 \ Bertiksichligt m man, daß Ai Abetangen von y, ae &, a Ya ao der angeführten Ab 
" lichkeitssatz fürGrenzschichten von derReynol ds schen Zahl unabhängig und beschränkt bleiben, 
so sieht man, daß durch die Veränderlichkeit von U mit x scheinbare Reibungsglieder in der 
 Störungsdifferentialgleichung hinzugekommen sind, wie schon 1929 bemerkt wurde. ‘Wenn mannur ; 
‚die höchsten Reibungsglieder berücksichtigt, entsprechend. dem asymptotischen Charakter der 


_ ganzen Stabilitätsrechnung für große Re ynoldssche Zahlen so erhält man in dimensiopslösen, EP 
/ Variablen, die mit U und 8 gebildet. wurden: KO BARS A 28 
k f | RR Te Fr 2: j E on 
Be „24v.. ay RU LEAD DD N 
er er dx ap Ray . 1% Be BER Be . Een 5 a u te NR: 3% 


Aka von U nach ® Ta: nicht mehr auf. Wenn man ‚den ee jeweils über BR 

x-Intervalle verfolgt, deren Größe mit b vergleichbar ‚bleibt, ‚zum mindesten weit kleiner als 7 

ist, das ja für die Änderungen von U mit x das Maß angibt, so kann man U als Funktion vny ie 
allein betrachten. Für wist dann ein Ansatz von a möglich, der auf die grund, 

legende 5 Drang een Bleiekliune, 1 Gl. 1 Male 


Alnepeansen‘ am 20.5.1944, fi % AH / e 4 Ban BR Er N E RR Be; 
_ Über die Diflerentalgleichung den transversalen. OR 
Stabschwingungen. | Be. 


_ Herrn Prof. Dr. L. Prandtl zum 70. Geburtsta 8. ee. 
Non F,P Jeiffer in Stuttgart. 


Es werden für die — die Querkraft und die Rotationsträgheit- berücksichtigende — ER partielle. 
Differentialgleichung 4. Ordnung der BREUER ale Stabschwingungen die Differentialgleichungen der 
Charakteristiken aufgestellt. 

Für die numerische Behandlung ı der Differentialgleichung erfolgt die Aufstellung der Ohemakieriptäbeh 
für ein mit ihr äquivalentes System.von zwei partiellen Differentialgleichungen 2. Ordnung und die Ein- 
führung der die Charakteristiken gebenden. Variabeln als neue unabhängige Veränderliche. er.“ 

Es wird die Fortpflanzung von Unstetigkeiten im Moment und in der Querkraft mit EIStjeN A 
der Charakteristikengleichungen beurteilt. Ba: 


Are set up the differential equations of the characteristics concerned imstehd of the linear Auetlat ER N 
ferential equation of the 4th order of the transversal beam vibrations taking into consideration ‘he tranwerse 1 
force and the rotational inertia. Fe 

With regard to the numerical method of the differential equation the setting up of the characteristics a 
will be made by means of a system equivalent to-the equation concerned of two particular partial differential Mn. 


equations of the 2nd order and by introducing the variables involving the characteristics as a new and ab- ED 
solute arbiütrary inconstant. 

T'he transmission of inconstancies in momentum and in the transverse force is gauged on the strength - 
of the characteristics equations. 


Sont posees les equations differentielles des caracteristigues en question au lieu de l’equation dif- Va 
ferentielle partielle et lineaire de la de ordre des vibrations transversales de bäton, prenant en considera- j 
tion la force transversale et l’inertie de rotation. 

Par rapport & la methode ‚numerique ‚de l’equation differentielle les caracteristiques sont posees au RK 
moyen d’un systeme equivalent a l’equation en question, c’est ü dire de deux Equations differentvelles par- Sa 
tielles de la 2eordre et en introduisant des valeurs variables involvant les caracteristiques comme variable 
nowvelle et independante. ER 

La propagation des inconstances au moment et dans la force transversale est r&alisee au moyen des 
equations de caracteristiques. 


Dun ımHeiinoro yacTHYHOorO ypasHeHna 4-ro paspana TOoIepe4HBIX IITAÖHBIX KOlTe- 
Öannf, IIPHHHMAIIIeTO BO BHUMAHHE NOIEPeUHYIO CHIY A BPAIMATENBHYO HHEPUNMD, YCTaHa- 
BIINBADTCH NuhpepeHmMHansHBIe YPABHeHuA XAPaKTepHCTHuR. 

Asa unctenHoro oßcyazeuun. nuhhepeHNMoHHOrO YPaBHeHHA IIPOHSBONUTCH YCTAHOB- 
_ HeHHe XAPäKTePHCTUK MIA ISKBHBAJIEHTHON C He CHCTEMH ABYX YACTHYHBIX AUPbepeHIHOHHBIX 
ypasHeHunü 2-TO paspana u BBeNeHHe IIeEPeMeHHEIX BEeJIMYNH, IIPHNAWIIMX XAPAKTePHCTHKaM 
HOBOe He3aBHCHMOe U3MEHEHHE. 
Pacıpocrpauenne (pa3MHo»keHue) TIepepEIBOB B MOMEeHTe U B NofepeuHoä cHIe YCTa- 
HABIUBAETCA IIPH IOMOIIH XApakTepHcTuK-BEIpaBHeHnd. 


Herr W. Flügge hat in dieser Zeitschrift, Bd. 22 (1942) in einem Aufsatz über ‚‚Die Ausbrei- 
tung von Biegungswellen in Stäben‘‘ die wohl zuerst von S. Timoshenkot) angegebene vollstän- 
dige a für die Transversalschwingungen eines Stabes 

022 02 0?k? 02 
2 i EI 
Ei te ekz (1 + Han AG 9y* 


ı) S: Timoshenko: Phil. Mag. Bd.41 (1921), S. 744—746; ebenda Bd.43 (1922), S. 125—131. 


(1) 


ne die Fortoflanzun N R\ 
PN Die Auflösung dieser Differentialgleichun ollt ufg 

> Indem vorliegenden Aufsatz sollen einige Hinweise gegebeı 
Einzeln diese! Differentialgleichung Aut Ans I Char | 


| “en. der able Ditferentialgleichung ist angenommen, AN si der 
i Er ströökt, ‚.yist die Zeit, z(x,y) die in der Biegungsebene senkrecht zur 
 biegung, e ist die Dichte des Stabes, q die Fläche des Querschnitts, | 
' Querschnitts für eine zur Biegungsebene senkrechte, durch den Schwerpu )) 
gehende ‚Achse, E der Elastizitätsmodul, @ der Schubmodul und A eine Kons ‚di 
Gestalt des Querschnitts abhängt. Später werden noch die QuerkraftQ und das Biegungsmom 
auf die sich die im Querschnitt angreifenden Kräfte reduzieren, und die Drehung ® de: 
„elements um die zur Biegungsebene senkrechte Achse durch den En 24 Qu 

Ne nz EB, G, gqund k seien lan BEN a 70 a 


1 Setzt ı man Hi So A z TE Meine ER 


une \ h AVzr , N £- L s £ Mr if 


x) PN 92 & 39 & vs i 02 a Br g 022 2 Wale BALIUEEN At RR Br 
A Be es re . EI ag PETE ER RE 
ER iR au. BEL Re: REN a A Dr 
Nr Va! FR! Br 7 oz N ; / 2 Ei = 02 N : = 2 = e 2 .n E . £“ 
EN ER N da u: en er ‚Om2dy2 ER a ray 37 rohr 
3 und | führt die Abkürzungen ein, RE AIPER NE Is EDN | a 
2 %2 er Eee 
2. k2 DE RT Fra, a TEE tr. 
h ER x 0 ai tag) En \ Ke, £ SH i 0 I EN j . % 1% Ka 
80: Bi die. De ee (1) © ae. ER RT h 
a | N en | nr FERNE NE SR TG RO 
ER Sie ist eine ae homogene Differentialgleichung 4. Ordnung und bildet mit den Gleichungen. FR 
iR B=mdetrndy, dß=rdetrdy, dy= de an udy, dö= det r,dy ö 
a ö — 
ein System von fünf linearen Gleichungen für die Größen a T, Tor T;. u der Deter- >. 
minanten der Matrix Se | ER BREI? 
a 0 —b 0 € Sy | 
da dy 0 DE EI0, 
HNO Yin de y 0 () n 
.e TR N) deze. dy 0 1 
0 0 0 de ..dy 2 


liefert für die Charakteristiken?) von Gl. (2) die Hauptgleichung > 


ady: —bdx?dy?—+ cdx!—=0 
und die Nebengleichung > | 


ade? (dade —dß dy) + (ady? —bda2) (dydc —dödy) + fdt=0 ...... (d).r> % 


| 
Aus Gl. (3) ergeben sich die 4 Werte für .- N 


%r = 'e, = VE = — @. 
f Y2 % E u er 
und mit der Bezeichnung x für einen dieser Werte wird Gl. (4): 
(da dß dy dö 
en A a ah Nie => K ! 
(5; #7.) + (az It 7.) +ft=0 RN RT 


Man hat also für dasjenige der 4 Charakteristikensysteme ‚,3. Ordnung‘, das zum Wert x gehört, 
die 8 Gleichungen: 


2) Siehe etwa R. CourantundD.H ilbert: Methoden der math. Physik, Bd.2 (Berlin 1937), 5. Kap. 


/ 


ein die IR Fokkonäh. yo: a BL ER he Y von Re N K% 
era Man kann! die 4 Systeme gewöhnlicher Differentialgleichungen nebeneinander angenähert 
“; „integrieren i in der Art, die Mass au?) für partielle Differentialgleichungen 2. Ordnung und Systeme ER 

' von partiellen Differentialgleichungen 1. Ordnung angegeben hat. Ich setze das Näherungsver- 


ER - fahren hier nicht auseinander; sein Grundgedanke tritt auch in Abschnitt 3 hervor, und die nume-, FR 


"0 Tische Durchführung wird dort ganz wesentlich einfacher. RS: 5 
A % Wir sen statt us: en leleietnug 1). das System der beiden Ditferentialgleichungen BR 
Ei Be E = 022 N 022 1) R N NEL Eee REN Mr ar 
RN Ne 4 a A | BIN te u 
Ba Dr Le N: Y ar SU, in | 5 g E Be, Jar En NE . : % a8 
TE TG Mr NT u na ae | ED SE ee 
zugrunde, aus denen durch Elimination von © die Gl. (1) folgt. ‚Sie En sich‘ aus den Gleichungen 
| FE B ar re. | x, = o Pe 2Q 22 £ BR ALERE 
ae et a gut 94% u Me 
F ee Kesujen od ERROR | 
SM gr ER und 9- _ 40 ee 
‚durch Elimination von M und Q*). Setzt man RT 
i ' N; 02 — x de 98 — ! 98 
E f . Bir y a7: 2» öy 9 M 
| | ee an OR NDR | 
ae Beau N aya oa a Buy ws 
‚so wird Gl. (6) RR ei PT ERS ee. 
\ ’ N ot, — /Gr, == — AGP, - } 1 (6°). 


ok?t, — Ek?r,— AGp, — AG® 
. Die Gl. (6) bilden mit den „»Direifenbedingungen ! er | 
dp, = r.de ze sıdy, dy = s,de +t,.dy, dp, = r,de + s2dy, Sa = sde+ttdy.. . . (7) 


ein A von sechs linearen Gleichungen für r,, 81, t1» "2 82 ba. Durch Nullsetzen der Deter-. 
minanten der Matrix 


re) 0 0 0 0 A, || Be" 
DREI OR OEL. 0 ok? AG (pı —®) a: 
de  dy 0 0 0 0 dp, 
0 de dy 0 0 0 dq, y 
Er 0 0 0 dx .dy 0 dp; \ 
I) a0 El! 0 da dy, dg; | 
ergeben sich die gewöhnlichen Differentialgleichungen für die Charakteristiken des Systems Gl. (6). 
Mit den Abkürzungen 
A@_,, VE=an net. RENTEN RIER TEENS 


| 3) J. Massau: Me&moire sur l’integration graphique des equations aux derivees partielles. Gand 1900, Chap. I 
u. II. Siehe auch Enceyklopädie der Math. Wiss., IIC2, Artikel Runge-Willers, $.161—164. 2 ER 
*) Siehe Handbuch der Physik (Springer), Ba. VI, S.359. Dort muß in der zweiten Gleichung Zeile 15 die linke 


0? 
Seite heißen: ogq%k? Ge" x 
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| werden est Differentialgleichungen für die 4 Systeme von Charakteristiken 


Hauptgleichung \ RER, Nebengleichung 
t: ee ; u. 0), 
1 e — Mr we. er - en +4AM a 0... (9b), 
ee er ee -.. | ip: A 492 440-000, oo), 
IV. | =: | BEN: ar 1 Eu, 4ao- 0... (9d). 


Die Projektionen der vier Charakteristikensysteme in die z,y-Ebene sind die Geradenscharen 


I. + c0y=/4, IE x —C1Y=Ag Ill. +63 y=A/s IN, e—oy—h .. (10). 
e - 1 +q ze drei Differential- 
gleichungen für die vier Funktionen y, 2,24 g, von x, und Analoges gilt. für die Differential- 
gleichungen II, ILL, IV. „Wegen des Auftretens von M, oder also p,, in Gl. (9a) und (9b) und von Q, 
oder also p,, in Gl. (9e) und (9d) bilden die Differentialgleichungen (9a—d) ein simultanes System. 
In der von Massau angegebenen\Art kann man diese Systeme nebeneinander näherungsweise 
integrieren. 

3. Es empfiehlt sich aber für die numerische Rechnung noch die folgende Vereinfachung. 

Wir führen A, und A, als neue unabhängige Veränderliche in die erste Differentialgleichung (6) 
ein. Damit wird diese Differentialgleichung zu 


Man hat in den Differentialgleichungen I (9a) zusammen mit 


022 1/08 29 
lt) EN ee at en Ba ALTE (11a) 
Weiter führen wir, und A, als neue VRanHane Veränderliche in die zweite Differentialgleichung (6) 
ein und erhalten 
02® Garslz 02 
04504; Akte A) FE Te NE (11b). 


Die Benützung der zwei Paare von unabhängigen Veränderlichen in den Gl. (11a) und (11b), 
nämlich A, und A, in Gl. (11a) und A, und A, in Gl. (11b) bewirkt, daß beide ‚Gleichungen die ein- 
fache Form, die nur die gemischte zweite Ableitung enthält, annehmen. Der Übergang von den 
Variabeln A,,A, auf A,,A, ist nach GI, (10) sehr einfach: 


Aıtk=i th N | (12) 
j 03 (A) —A,) = 0, (Ay —A,) 
Führt man ein 
f od 1 
M=EgaW 77 Pat 9) R 
und ea ME (13), 
02 
= A64(5,—®) = A@4 (pı + 1 —P®) 
so werden die Gl. (11a) und (11b): | 
022 92® 
59,07, — AM und EIREN = A ET TAN (13a), (13b). 


Wir behandeln die Differentialgleichung (13a) wieder mit dem Charakteristikenverfahren. Setzt 
man 


d2 Dr 022 022 0?2 
a La, du) oe 
so wird Gl. (13a) 
= AM. 


und bildet mit 
dp, = rd, +3, d}, und dg, = s,d}, +t,dA, 


i drei linearen Gleichungen für 1 ne tu 


Be Are d r Matrix nt i 

naeh; N I U 0 h Ina RE 
WR 2 a a Da ALIEN I dA a, | 0° 1 .de, ET ar 
aan ld, Teak ınde 1er  ae 
5 Hl. b> gerinnt man für die Charakteristiken Sy U ER re NE ö 
RER i 2 der Schar TR #32 den Schar ILL" k 
edle nee Nora Ka SEN ‚8 ET SA CODEE TAT 
Ber wüdie, Nebengleichung . . . : dp, = AM di, (da, : dy=AMdA, | (14b) fer. 

I Wozu die Streifenbedingung tritt = 4 Ah rn | dz = p, dA, ö + 
Venen re Ganz analog erhält man zu Gl. (13b), wenn man setzt x 
De} od a ee a Le 
ee 9 er TEA ® Dr Ai = 
| , . für die Charakteristiken & : a A 

a ve ’ der Schar III der Schar V AN : 

die Hauptsleichung 4... 2). A, const, RL A, = const. Rat | Be. 
die Nebengleichung .. .. .d,y =-—AQdA,\ (14c) dg= —AQdi,! (14d) a 
wozu die Streifenbedingung tritt dd = 92 dA, d® =p,d}, Bi 
Für jede Schar von Charakteristiken hat man in ihren beiden letzten Gleichungen zwei Gleichungen = S 2 
für drei unbekannte Funktionen, z. B. bei ur Schar I zwei CHSEOHWERSR für 2, 2, Qı als Funktionen re 

yon 2. a 
K Weiter ist noch zu beachten, daß in en Giernngen für die Scharen Iund IIdas M,. oder 


also 9,+g, [wegen Gl. (13)], auftritt, das durch Integration der Differentialgleichungen (14c) Br. 
und (14d) gewonnen werden muß, und daß in den Gleichungen für die Scharen III und IV das @, 5 
oder also p9,+-9,, auftritt, das durch Integration der Differentialgleichungen (14a) und (4b) 
gewonnen werden muß. Zur angenäherten Integration der Charakteristikengleichungen in der R; 
Art des Massauschen Verfahrens muß man also die den vier Systemen [Gl. (lMa—d)] entspre- 

chenden Differenzengleichungen nebeneinander heranziehen. 


L 


Bild 1. a: 


Seien etwa die Werte von z, 9,91; ®, 3, 9, in den Punkten der x-Achse gegeben. Dann Be 
geht man (Bild 1) vom Punkt 1 der x-Achse auf der Charakteristik I (A, = const.) der Schar I Ber 


und von dem (hinreichend nahe bei 1 gelegenen) Punkt 2 der x-Achse auf der Charakteristik II 
(A,—=const. ) der Schar II weiter bis zum Schnittpunkt 3 dieser beiden Charakteristiken und be- 
rechnet p}, q}, 2? (die oberen Indizes geben, da Verwechslungen mit Potenzen ausgeschlossen sind, 
die Nummer des Punktes) in Punkt 3 mit Hilfe der den Gl. (14a) und (14b) entsprechenden 
Differenzengleichungen | 
} pi = AM a E—G-AMM—A)N nach | 
Ann r en e nach (Gl. (14a), 3 en nach Gl. (14b). 
Daraus erhält man p° und g?, ferner doppelt den Wert von 2°, wobei die beiden Werte 2? innerhalb 
- der geforderten Genauigkeitsschranken übereinstimmen müssen, sonst muß 2 näher an 1 gewählt 
werden. 

Durch 3 zieht man nun die beiden Charakteristiken III (A3=const.) der Schar III und IV 
(A,=const.) der Schar IV, die die Punkte 4 und 5 auf der x-Achse geben. Von diesen aus rechnet | 
man auf den Charakteristiken III und IV nach 3 mittels der den Gl. (14c) und (14d) N, EN 

chenden Differenzengleichungen 


—p = —AQ*(R) | ,—g= A—N 

3 Ne r nr nn 4) nach Gl. (l4c), g= apa an ») nach Gl. (14d), 
Daraus erhält man p} und g; und doppelt den Wert ®3,. Hinreichende Übereinstimmung der beiden 
Werte ®3 muß allenfalls auch durch Verringerung des Abstandes 1—2 erreicht werden. 


> 
1 
% r 


ei Luntikige ae 


und ABA = neue Werte a () und Cs 


\ 


Von 2 und 6 aus ( Strecke 2—6 gleich Strecke Ss wird Punkt 7 gefunden als Schnitt A 
durch 2 und 6 gehenden Charakteristiken Iund II. Von 7 aus werden 8 und 9 mittels der Charak- 


: De y>“ 
1 ai . 
y P 


teristiken III und IV durch 7 bestimmt und von 2,6, 8,9 aus die Zustandsgrößen in 7 berechnet. 


In dieser Weise ergeben sich weiter die Zustandsgrößen in den Punkten 10 und 11 in ‚gleicher H 
Höhe über der x-Achse wie 3und 7, und mit diesen Punkten ist das gestrichelte Netz mit seinen 
Schnittpunkten auf III- und IV-Charakteristiken bestimmt. Die Zustandsgrößen in jedem Punkt 


der Strecken 34, 3—1, 3—2, 3—5, 7—2,... . werden dann als bekannt (durch lineare Interpolation 


Charakteristiken III, I, IL IV von 3, 13, 14, 7 aus usw. 


In dieser Weise lassen sich die Werte 2, P,, qı ©, Ps, 43 berechnen in diskreten Punkten 
Punkten eines Dreiecks in der &y-Ebene, das begrenzt ist von einem Stück AB der x-Achse und 


den beiden Charakteristiken III und IV durch A und B, wenn wir voraussetzen, daß III u IV 


gegen die x-Achse flacher verlaufen als I und II. 


4. Wir geben den Gang des Rechenverfahrens noch en den la Spezialfall von Rand- 
und Anfangsbedingungen. 


Der nn erstrecke sich von x=0 bis =. Zur Zeit y=0 sei 21 0) = 0. Damit ist 


für y=0 auch = identisch in x. Ferner sei 0 für an und alle x, ®(x,0) =0 und damit 


od oDd 
auch 0 für y=0'und alle x. Schließlich sei 0 für y=Oundalle x. Alles entspricht dem 


EHnstand, daß der Stab zur Zeit y=0 im Ruhezustand sei. Damit sind auch M und Q zur Zeit 
y=0 für alle «= 0 gleich Null. In dem Gebiet rechts der Charakteristik IV durch 5 (Bild 3) ist 


alles Null. Das entspricht der Anfangsbedingung und den Differentialgleichungen (13a) und (13b). 


' zwischen den Endpunkten der Strecke angesehen. Die Zustandsgrößen in 12 ergeben . auf x h 


Längs x—=0 sei 2(0,y)=0 (Stab bei «—=0 dauernd festgehalten) ; damit ist 0 längs 


‚oD Ä ER 
der y-Achse. Ferner sei 3,0 y)=/(y) eine gegebene Funktion von y mit f(0)=0. Es wirkt also 


von y—=0 ab ein im allgemeinen zeitlich veränderliches Moment an der Stelle «=0 auf den Stab. 


ö In 6 (Bild 3) sind bekannt: z=0 und am). Also gilt in 6: 2°=0; pi=q!. 
Wegen ap (y) ist, weil — 3 =; AN p% + g° = f(y®). p° und Dserhält man auf III von Er 
aus mit on den Gl. (14c) entsprechenden Differenzengleichungen‘ 

9%, —m, = —-AQW—A) 

Deniz): 
Damit hat man auch g} aus 95 = f(y®) —p}. 


N erhält man ; ie ala der zweiten. or (14a) I Ditferenzengleichung A 
er : M—PM = AM” (4% — A!) längs I von 7 aus, wo die Zustandsgrößen durch lineare Interpolatieuue 
zwischen 5 und 3 bekannt sind. Damit hat man nun z, Pi» 4 PD, 92, 9; in 2 RL 
MR ur Die Zustandsgrößen in 8 erhält man auf: N BR 
Br - Charakteristiken IV, II, I, III von 6, 9, 10, 11 VS FR Ten 
aus, wobei die Werte in 9 durch lineare Inter- | 
..- polation zwischen 6und 3 erhalten sind und a 
der Charakteristik IV von 5 über 7, 3, 10, 11,. 
alles von rechts her bekannt ist. 
42 Damit kann nun der Streifen zwischen 
den Charakteristiken IV von 5 aus und IV von 
6 aus ‚erledigt ‚werden. .Dann wird Punkt 12 
' analog wie 6 behandelt und der Streifen zwi- 
schen den Charakteristiken IV von 6 aus wu 
IV von 12 aus erledigt usw. 
00.5... Der Nachteil des, Verfahrens ist die un- 
günstige Fehlerfortpflanzung. Sie bewirkt, daß 
R - man trotz Verkleinerung der mit 1—2 bezeich- 
' neten Strecke und der damit verbundenen Ver- 
mehrung der Rechenarbeit bei weiterem Fort- 
schreiten in der y-Richtung immer wachsende 
Unsicherheit des Resultats bekommt. ET 
0.8. Wie weit gibt das Näherungsverfahren 
-  Aufschluß über die Fortpflanzung von Unstetig- 
.. . keiten im Moment M und der Querkraft Q? 
= Es ist 


w 


Fe | 
EN 


M-Ege — Ban, 5, 
0% BR 
02 a ZEN 
. Wir setzen voraus, daß M nur im Punkt 1 einen Sprung (AM)! habe. Damit hat p, in 1 einen i A 
Sprung (A9,)!= ne (AM). Alle andern Größen seien auch in 1 stetig vorgegeben, insbesondere 2 
‚ist also AM). : | en 
Wir gehen nun auf die Variabeln 4, über. Nun gilt 
PR=P+9 und %=0(p—g): 
Alsoist | MAMA) > 
und wegen (A9)!=0: “ Ber 
| (Ap2)t = (A)! 


als Folge unserer Voraussetzung, daß nur M bei 1 den Sprung (AM)! hat. 


1 ..: 
Wegen M= ZH (P2 + 9:) 
; | 
i RE TEE 1 5 
: Längs der Charakteristik III von 1 aus gilt nach Gl. (14c), wenn wir die Differentialgleichung 
eh eine Differenzengleichung ersetzen (Bild 4) NS 
m = ARM —R), SE 
und zwar beim Fortschreiten am linken Ufer iz’ DS 
von III FR ET 3 r 
A \ 7 74 
Pa — Pa = — ANA, —A,) £ BIS 
NR END N UENERN DE 
und am rechten Ufer - nz x » 
u ee HE ———— 4 — 
? 2 2 f Val [9 4 


Hieraus folgt 
(An)? = (Ap)!=2A(AMN)!. Bild 4. 


3 x 


90° Pfeiffer, Über die Differentialgleichung der transversalen Stabschwingungen ee 
Für q, gilt längs der Charakteristik IV von 2 aus nach (14d): 

= AP, —A), | 
also, wenn wir wieder am linken und am rechten Ufer der Charakteristik TV von 2 nach 3 gehen 

1 —B= AR (A, —A,) 

| a —n= TARA). 
Hieraus folgt 
(Ag’=0. 
Für (AM)? erhalten wir also 


1 1 
(AMP = AP)’ + (4) = (AM. 
Wir beachten noch, daß Q? keinen Sprung hat, denn p} und g; werden auf Iund II von zwei Punkten 


her gerechnet, wo keine Sprünge sind. 
Wir rechnen noch einen. Schritt weiter auf III von 3 nach 5. 


Nach Gl. (14e) ist pP —p} = —AQ(A} —A,), 

und zwar Pr — or = AP (N — A) 
und AR). 
Daher ist * (A9)°=(4p)?—=2A(AM). , 


Für 9, gilt längs der Charakteristik IV von 4her _ 
B—g= Ag —h), 
also wieder h (Ag. 
Q5 hat keinen Sprung, denn es wird auf Charakteristiken, auf denen keine Sprünge zwischen 


linkem und rechtem Ufer sind, von anderen Punkten her gerechnet. 
Für (AM)? erhalten wir 


1 1 
(AM)’= 1A (42,)’+ (Ag2)]= Dx (AM). 
Es ergibt sich also (beim weiteren Fortschreiten): Der Sprung (AM)! pflanzt sich zur Hälfte in 
konstanter Größe . (AM)! längs der Charakteristik III von 1 aus fort. 


Die ganze analoge Rechnung zeigt, daß sich die andere Hälfte des Sprungs (AM)! in kon- 
stanter Größe längs der Charakteristik IV von 1 aus fortpflanzt. 
Setzen wir nun voraus, daß nur in 1 ein Sprung (AQ)! der,Querkraft Q@ sei. Damit hat 7, 


& 1 e : s x 
einen Sprung Adg (AQ)!. Alle andern Größen seien auch in 1 stetig vorgegeben, insbesondere 


ist (Ag,)'= 0. Nun gilt 
PA=-9+g und 1=0(Pı —4)» 
also IR Aq= c,(Apı —Agı). 
Daher ist wegen (Ag,)!=0 wieder 
(Ap)'= (Ag)! 


als Folge unserer Voraussetzung, daß nur @ bei 1 den Sprung (AQ)! hat. Also ist 
(AQ)"= 2 A@g (Ap)'= 2 AGg (Aqı)". 


Eine ganz entsprechende Rechnung zu der, die die Teilung und Fortpflanzung eines Sprunges 
von M zeigte, ergibt, daßsich ein Sprung in Q an der Stelle 1 zur Hälfte längs der von 1 ausgehenden 
Charakteristik I und zur Hälfte längs der von 1 ausgehenden Charakteristik II in konstanter 
Größe fortpflanzt. 

Wir behandeln noch den Fall, daß in 1 das Moment M einen Sprung (AM)! und die Quer- 
kraft.Q einen Sprung (AQ)" hat. Alle anderen Größen seien auch in 1 stetig vorgegeben. -Damit 
hat 9, in 1 einen Sprung (Ap,)! und 7, in 1 einen Sprung (A7,)!, aber es ist (A9,)!= (Aq,)!= 0 
vorgegeben. 


Wegen na n=4(m—q) 
gibt (Aq,)!=0: (Ap)'=(Aqı)l 
und wegen 92 6, (Pa —9;) 


gibt (Ag,)!= 0: (Ap)'= (Ag,)!. 


A 


y ne] Yuan OGK n ar 
amı= 2% (Andi dag 


ET d ED LrgaS (OR 2 AR (Ami Ahren), 0 Aa 
FRE N Rechnen wir jetzt längs der Charakteristik III von 1 nach 3 in a, so gilt Wegen Gl. „di 40): #:: 
BURN. Dee AN a 
A R us = a y | Y h Eu 2 vr, Par — Par = = Ag, ( h ERS: 8 Pax Äi IR AR 17 x TE ö 1 nd 
Re a RR pa) (AQN. P% h 


2 ee Fi Und rechnen - wir Ines der Charakteristik IV von 2 nach 3, so, ein Gl. a 
$ ö  mM—R=—ARM—R) 5 5% 


ie 


5 = “ 5 2 2 Br aa er = —4Q2 (2 —A?), A N 2 
ee EL ee A De a | 
Daher u Er An Re | HR 5 

. (Amy=,llApJ) + (49) Aa A91= Lam ı ( 2) ao. 


(A rOrR 2 p und. q, längs ek I und II Yon! a her gerechnet werden, an 
nn N für alle Größen vorausgesetzt ist. Damit ergibt sich nun für u 5 ie Ib, 
„0 Nach Gl. (14c): a nach Gl. (14d): 
SEHR a Er A Te 
en BUN 
(dq)’= a i | 


= Aa A) 
u 
(Ap,)°= (Ap;)? 


Aa A 5 ” 


Also | Ka 

am= And + (da) (Apr al 12) (40)! uw. 
Da A —A, beliebig klein angenommen ae “en so zeigt die le daß Se die Hälfte . 
des Sprunges von Minlin konstanter Größe N m von 1 aus längs, der Charakteristik III 


fortpflanzt. Die analoge Rechnung zeigt, daß En die andere Hälfte von (a M)! von 1 aus Tanne 
der Charakteristik IV in konstanter Größe fortpflanzt. 
| Rechnen wir von 1 auf der Charakteristik I nach Gl. (14a)] und von 1 auf der Charakie- 
ristik II fnach GI. (14b)], so erhalten wir in ganz analoger Weise, daß sich der Sprung (AQ)! von Q 
in 1 je zur Hälfte auf diesen beiden Charakteristiken fortpflanzt. 
Von einem Sprung (A M)!des Moments M und (A Q)!der Beer Q in einem Punkt pflanzt 


1 
sich also ie — (AM)! längs der Charakteristiken III und IV und je - 54 une der Charakte- 


ristiken I En II fort. 
Eingegangen am 2.Febr. 1945. 


KLEINE MITTEILUNGEN 


Die Wirkung der Puffermasse beim Zusammen- 
stoß von Eisenbahnwagen. 


Stoßen zwei mit: gefederten Puffern versehene Eisen- 
bahnwagen zusammen, ist die Schwingungsbewegung 
leicht zu berechnen, wenn die, Federn keine Masse 
haben [1}). Die Masse der Puffer und Federn ist klein 
gegenüber der Wagenmasse. . Bei einem, von der Reichs- 
bahn viel verwendeten Puffer ([#], S. 992) haben die 
bewegten Teile (Pufferhülse mit Teller und Stoß- 
stange voll, Pufferfeder zur Hälfte gerechnet) ein 


Gewicht von etwa 92 kg, beineueren Pufferbauarten [2] 
nur etwa 6ökg. Beim Zusammenstoß von zwei Wagen‘ 


arbeiten vier Puffer gleichzeitig, die bewegliche Puf- 


4.92 
fermasse ist also z. B. iz 0,375 kgs?/em, die Masse 


eines Wagens dagegen 25 bis 100mal größer. 


Gottschalk [3] hat beobachtet, daß beim Auf- 
laufen eines schweren Wagens auf einen leichten die 


1) Die Zahlen in eckigen Klammern verweisen auf das 
Schrifttumsverzeichnis am Schluß der Arbeit. 


Puffer des auflaufenden Wagens bei allen Auflaufge- _ 


schwindigkeiten sich mehr zusammendrückten als die 
des Prellbockwagens. Wir erklären diese Erschei- 
nung dadurch, daß wir die Puffermasse bei an Zu- 
sammenstoß berücksichtigen. 


1. Beliebige Massenverhältnisse. 


Wir nehmen an, der Wagen a mit der Masse M — Ma 
fährt mit der Geschwindigkeit v, auf einen stehenden 
Wagen c mit der Masse Me= M :n auf. Zwischen 
den Wagen befindet sich die Federmasse M— = M: m, 
zwischen M, und M» und zwischen My und M, je eine 
Feder mit der Richtkraft c. Die Zeiten zählen von der 
ersten Pufferberührung ab, die Wege von der Lage 
der Massen zu dieser Zeit (Bild 1). Wir betrachten die 
Bewegung von einem Koordinatensystem aus, das sich 
mit der Geschwindigkeit des Massenschwerpunktes 0) 
bewegt: 


OR mn a N Wer 
I mntm+n (1) 


ET Bi 1. Wederschaltbild. 


"Nach don ee Er Anfangsmustand. erchten An 
Sc nahme setzen wiran: 


Be — (4 sin wi,” 
De = Ce sinwi 


j 


= Obs sin 07 


zung durch sin wt: 
n. cC 


20) Op+e0, =0. 


uw —e) 0 Er 
'/M 
m 


| S 


en 


+ wo 


der ne W=Y4m+(n — 1%: 


Wr VairyamFa Tr 


2 
Br » r £ 18° 228° 
> . } 
Be." 2. Bild 2. 
% En . Das Verhältnis 
ER R DE N . (6) 
ist nur von m und n abhängig. | 
i %a = Oaısinwt+ OasSinwst. (7) 
e Die Konstanten C, und Oo rechnen wir mit den 
Gl. (4) aus. 
Ca a1 s i 
B»=— 9 (2?m+n—1— W)sinwt \ 8) 


1 


> (2m+n—1+ W)sinw,t | 


NC, 
77 


"(na —1— W)sinw;t | (9) 


Car 


+ m—1+ Wsinwt | 


Für t = 0 ist die Geschwindigkeit der ‚Wagen. unter 
Beachtung von Gl. (1): 


m+n 
Bug Yo mn tmtn’ . (10) 
d Z MW 
id = den = — Yon mEn 


R Bo weile P abmax, und an 
den Größtkräften. Bes die beim Vernac 
Puffermasse Bu Sa les N “ an 


SEN 
SRSREL nebeneinander und. hintereinander ‚ges d, 
während ce für zwei a. Fedeı ar 


Ei t al. 2 und nach Ken) 
und finden durch Einse zen. in. ( ) | h eilt: so daß wende ist: 


a 


Die Hauptdeterminante des Ausdruckes: (4) muß ver- 
‘schwinden, wenn eine endliche Lösung existieren soll. _ 
Nach Kürzung durch die triviale Lösung w® = 0 ist mit 


Sn 


cos W=+ 089 u. ; 3 IS 


in “4 


Ze 


an! 


f = 2, ER 
Ss % en u 


PB it Hi Ba re a von a 
vier Pufferfe ae and 2% 


x 5 : De =% Veen: = Fin. BE 


= re 


Er der g- Pabmax: Pax Fr RE) 
ande Ne N Ne 
ER ER me Phomax: Pranke ae 14) > 
2. RK ANe vr L RAALE 


Wir nehmen M„=30, My» 0,395 und M;=. is 
2 Bere an (m = 80 undn ='3). Nach GI. (5) wird. 
—1,412 YeIM und w, 12,728 a und nach“ 


Gl.(6) p= 9,014. Nach Gl. (11) wird 0x, w, =, 26203 © 
und Ca, w, — 0,00314 vg. Die Gleichung ' u = 0 En 
lautet, wenn man w,t = u setzt, : 


008 1 0.9750 con D.OIA a0, 


Die Gleichung hat verschiedene Null-Stellen. Der N 
größte Wert für % ist in der Nähe vonu=n2? um 
erwarten (Bild 2), und man findet die Lösung für 

u = 89,88%. Damit wird P45 max = 0,4100 v, at 

Die Gleichung te = 0 lautet: 


c08 u — 0,9997 cos 9,014 u — 0. u 


Die Gleichung gibt bei u — 89,88% für &.e einen 
Kleinstwert. Die benachbarten Höchstwerte liegen 
bei w, = 71,90° und u, = 107/85°. Die zugehörigen 
größten Kräfte sind Py. = 0,3859 v, VelM und 
0,3864 v, Yc/M. Die Feder a b wird also um 6,1% 
mehr beansprucht als die Feder dc. Nach Gl. (12) ist 
die größte Kraft ohne Berücksichtigung der Puffer- 
massen Pihax = 0,3536 vo YelM, q=1,160 und 
n—= 1,093, dh. die Federn werden tatsächlich 16,0 


bzw. 9 3% mehr beansprucht, als sich bei Vernach- 
lässigung der Federmasse ergibt. 


2 


ieh 


0 8.DarstellungderV 
ae ua iue nel, . R FIRE 
"Haben der stoßende und gestoßene Wagen gleiche 
ER Massen (n — 1), vereinfachen sich die Gleichungen. 


erhältnisz 


Fr H p=YV2m+1 Pr Be ERNST NH RS He FE 


In. M : 2 er ER 
Os = 3% V, ud 0,=0s,:7° (Me) 


Die Amplituden der Oberschwingung sind also im Ver- . 


... ‚hältnis 1:9p® kleiner als die der Grundschwingung. 
52 Die Pufferkräfte sind  . 
A, ar $ F { 
Po $.3 RER j 1 Nr 
Par -=7 V:M (sinu+ 2 sin pu), 
x N nd do. ASIEN eu \ 
BR ET Pro=zyYeM (sinu—sinpu), 


701 70 ER] 


0 0 0 MM WO 120 
; Bild 3. Verhältnis der größten Pufferkräfte. 


- Den Höchstwert für P,» finden wir aus der Bedingung 


cos u+ cospu=0. w=(22z—1)n—pu, wobei 
z=1,2,3,...undzmax + + 1 ist. 

@e—)a 

ERS 
- Dabei ist sin u = sin p u und damit 
r 
1\ı » 2@2—]1)r 
a (147, )sin D+l re (13a) 


Für den Höchstwert von Pp. gilt die Bedingung 


cos u — cos pu = (0, ER wobei 2 =], 
RT. 1 nr N: 
2,3,... und zmax SEN U 91 . Dabei ist 
sinu =-—-sinpuund 
DElav 9 2762 
r=(1+7) ET C (14a) 


Die Zahlen g und r sind in Bild 3 über p dargestellt. 


Di x m-Vr und wW— mVm+Ll. Ba) 


Re % ist die Kreisfrequenz der Schwingung aka Puffer- 
masse und w, die Kreisfrequenz einer Oberschwingung. 
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EINÄT, i 


DEREN RN ER 
' Die Hyperbel 14 hüllt die g- und r-Kurven ein. 


Die größten Werte von q und r sind immer größer als 


‚1, d.h. die Berücksichtigung der Puffermasse erhöht 
. Immer die Federkräfte. Im Bereich der im Eisenbahn- 
‚wesen vorkommenden Verhältnisse ist m etwa =25 
‚bis 160, d.h.p = 7 bis 18. Die Beobachtung Gott- 


schalks, daß die Feder des auflaufenden Wagens 


mehr beansprucht wird [3], bestätigt sich in denBe- 
reichen p = 8 bis 10, 12 bis 14 und 16 bis 18. Inden 


Bereichen p = 6 bis 8, 10 bis 12 und 14 bis 16 ist das. 


> 


‚Gegenteil zu erwarten, hier wird sich die Feder des Ir 


angestoßenen Wagens mehr zusammendrücken. 

Zeh >, ur > ee ! 
4.Ergebnis. Re Marsch DENE 
Vernachlässigt man beim Berechnen des Zusammen- 

stoßes von zwei Eisenbahnwagen die Masse der Puffer- 
federn, erhält man größte Pufferkräfte, die bis 15% 


: kleiner als die wirklichen sein können. Die Puffer 


des -auflaufenden und des  angestoßenen Wagens 
drücken sich im 
sammen. 


Schrifttumverzeichnis. 
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(1934), H.23, 8.421. 
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Bd. III, 26. Aufl. Berlin 1934. 
Oppeln. GerhartPotthoff.'. 
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Über das Abklingen von Schwingun- 
gen mit schwacher in beliebiger Weise 
von der Geschwindigkeit abhängiger 
Dämpfung. 


Wir betrachten eine Schwingung, von der 
Form 


st /@)tor=0 ....(d) 


wobei f(x) in dem durch die Anfangsampli- 


tude a, festgelegten Intervall— a, St <S 
+ ®a,eine ungerade analytische nicht ab- 
nehmende Funktion sein soll. Wegen der 
Approximierbarkeit allgemeinerer, etwa ge- 
knickter Kurven durch Polynome ist dies 
keine wesentliche Einschränkung. 


Der Betrag von f (x) sei so klein, daß während einer 
Schwingung die Amplitude nicht merklich abnimmt, 


also | f (&) | << w?:ay. i | 
Multiplikation, von (1) mit & und Integration über 
das Zeitintervall zwischen zwei konsekutiven Null- 
stellen von % mit den zugehörigen x-Koordinaten 
= — 1, und = + 2, ergibt 
w* 


o+[taadı+ 7 e—eh=. 


Gemäß der vorausgesetzten geringen Dämpfung und 
der hieraus folgenden langsamen Amplitudenänderung 
führen wir die Amplitude a zu Beginn der einzelnen 
Schwingung und den Betrag Aa der Amplituden- 
änderung während einer Vollschwingung (da >o) 
ein, so daß unter Vernachlässigung von Gliedern 


höherer Ordnung geschrieben werden kann 2, = — a, 
Aa ® 
% = a— —— und 
2 


allgemeinen verschieden stark zu- 


rtechnischen Mechanik, 


94 Kleine Mitteilungen 


1 n=—aA a, 
[t@ ad — o:a4a=0. 


Weiter kann wegen der schwachen Dämpfung das hier 
auftretende Integral genähert aus der ungedämpften 
Schwingung berechnet werden, es darf also 


x = wa? — a2 
gesetzt und als Integrationsintervall die halbe Schwin- 
gungsdauer 2 T, genommen werden. Es ist dann 


a 
ada=L . free = 
0) 


RER: 
2 


a 
| 

0:10 = | tolar—as) dr 
—q 


a 
> 
ham | 1aYaaayda. 
0 


Ist die Auffassung der Amplitude a als einer stetig- 
differenzierbaren Funktion der Zeit gestattet (Hüll- 
kurve nach Ze chl), so kann weiter geschrieben werden 


lat 
Asa nu, 


R © 
RG. 
also 


a 

3 2 RE REN. 

a-— | tele ande (2) 
0 ! 


oder mit der Substitution & = a:cos p 
7 


2 

. 2 

oa=— — | fo asin p) sin pay 
0 

Die explizite Berechnung des Amplitudenverlaufs für 
gegebenes f nach (3) läuft jetzt ersichtlich auf die Be- 
rechnung eines bestimmten Integrals hinaus. | Wir 
wollen hier jedoch die umgekehrte Aufgabe behandeln, 
aus dem beobachteten Abklingvorgang auf das Dämp- 
fungsgesetz zu schließen. 

Es sei also die Funktion a =«(f) und damit auch 
a(t) gegeben. Man verfügt somit über eine Parameter- 
darstellung der Funktion a (a), die wir im folgenden 


benötigen. Die Funktion f(x) soll nun berechnet 
werden. 
Hierzu führen wir auf der rechten Seite von (2) die 


neue Integrationsvariable u ein: 
u=ot@ a), a=]/a—M, 
’ fer} 
d 1 ..du 
TR SE 
a — — 
2 


und erhalten 


Ba) 

abo 

a=— 7 an 
ro a) Yora—u 
oder mit v— w?a? anstelle von a 
® 
2. (Almen a 
Re 7L u a ei ee ) 
0 


2) Th. Zech,Ing.-Arch. Bd. 13 (1942), 8. 21—33. 
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Dies ist eine A belsche Integralgleichung für f. Sie 
hat die Lösung?) 


U 
u v 
0 

Schreiben wir jetzt x anstelle von / u und ersetzen wir 
die Integrationsvariable v wieder durch a, so erhalten 
wir wegen v— 2 w?aa die gesuchte Dämpfungsfunk- 
tion f (&) in der Gestalt . \; >> 

zo. 

d a(a)a?da 


a 0720, nor DE 
F 


Besteht speziell zwischen a und a eine Beziehung der 


Form a = — c-a”, so tritt in (5) auf der rechten Seite 
ein Integral der folgenden Form auf: 
1 1 
ırdx == et rd 
——— I (1-22) ?2rdr=3 |1--y). ?y\erd 
{) 


n+]1 

Iz(_in-l Kar =) 

Tag 22 TETLENS, 
3 


und es ergibt sich für die Dämpfung 


n-+ 3 
1 ar 2 r( 2 ) ee 


Dresden. Plato. 


Eingrenzen der Wurzeln von Gleichungen 3. 
bis 6. Grades. 


Wenn man sich rasch einen Überblick verschaffen 
will, wo die reellen Wurzeln einer Gleichung 3. bis 
6. Grades etwa liegen, hilft oft ein Verfahren, das 
hier am Beispiel der Gleichung 6. Grades gezeigt wird. 
Das Verfahren gibt mit voller Strenge darüber Aus- 


kunft, wo reelle Wurzeln ausgeschlossen sind; bei 


Gleichungen 3. bis 5. Grades kann man den zu unter- 
suchenden Bereich beliebig vorschreiben. 
Die Gleichung 6. Grades mit konstanten reellen 
Koeffizienten 
"rad ++ +dntes+f=0. (la) 
werde umgeformt zum Beispiel in 
(0° +c)(&® ax +b) >. (Ib). 
= (ac —d)@® +(be—e)x—f h 
Nur für die folgende Erklärung werde dafür kurz ge- 
schrieben 
1) OK) Fi Fear) PO Br (le), . 
worin e 
Y=3 »=Rte y=a+as+b (Mbis(d) 
und 9 = (ac —d) #2 + (bc—e) 8 — f gesetzt sind. 
Die leicht aufzulösenden Gleichungen 


y-9% »=08 y=-09, =0 
ergeben diejenigen reellen x, bei denen eine der Funk- 
tionen y das Vorzeichen wechselt, wenn x alle reellen 
Werte von — oo bis + ® durchläuft. Diese „‚Grenz- 
werte‘ seien der Größe nach geordnet 


HKHII rn (6). 
In dem Bereich zwischen zwei aufeinanderfolgenden 
Grenzwerten kann Gl.(lb) durch einen reellen Wert 
von & nur erfüllt werden, wenn die Vorzeichen von %, 
bis y, zueinander passen. In den Bereichen, wo sie 


?) G. Hoheisel, Integralgleichungen, 8.117. Sammlung 
Göschen. Berlin u. Leipzig 1986. : 
‚6. Kowalewski, Integralgleichungen, 8. 9. Berlin u. 
Leipzig 1930. 


1en, a tler uns 


eschlossen. Obin. 
turzeln danach in Frage 


ap sich ei Wurzeln befinden und wie 


eh +40 0° + 4410 720-0 ı 


A viele, muß dann durch a es Zahlenwerte 
' „untersucht werden. h 


Beispiell:, $ ER 
4305 —4l as 3% ie (7a) 


wird umgeformt i in 


(8-81) —4) Ir 
anal x + 3123 2 + 720 \ 2b) 
See ee SE 

REN ENT RNCNH 
Y. = ‚gibt 2 — 87 = 0 mit der reellen ' (0), 
Lösung s=+ 4,4 
a =0 bedeutet a ee) 1 . (10), 
und gibt s=+5luds=-—81 ae 
Ya 0 bedeutet # — 475 —11=0 E 1). 
und gibt = +49 und = —0,3 )J 
‘ Die Grenzwerte lauten geordnet 
ER 7 2.0.35 0 AU RAO ,T . (12). 


Für jeden Wert von & unter — 8,1,z.B. A =. 710, 
werden y, negatiy, Y, negativ, Y, positiv, %, negativ, 
somit kann hier keine reelle Wurzel der Gl. (7) liegen, 
denn schon die Vorzeichen passen nicht zueinander. 
Für jeden Wert von & zwischen — 8,1 und — 0,3 
werden alle vier y negativ, folglich sind hier reelle 
‚Wurzeln nicht ausgeschlossen. In Fortsetzung dieser 


Betrachtungen findet, man, daß in den gestrichenen 


Bereichen reelle Wurzeln ausgeschlossen sind: 
E82 208..20,.440 49 


Der schmale Bereich 4,9 <x<5,1l macht es leicht, 
durch probeweises Einsetzen © — 5 ‚0 als Wurzel der 
Gl.(7) zu erkennen. Bekanntlich ist dann minde- 
stens noch eine reelle Wurzel vorhanden, doch ehe wir 
versuchen, sie durch probeweises Einsetzen zu finden, 


können wir die Bereiche noch weiter einengen. Gl. (7 } 


läßt sich noch auf eine 2. Weise umformen: 


(2er 3 En, 70) 
= — 52 0° — 4445 + 720 
nach der Formel 
. (74). 


(0° +0) (Has Hb—cla))-, , 
=[(b — c/a) (cJa) —d] 0 —ex — f 
In Gl. (l11ec) schlagen, wenn % alle reellen Werte durch- 


läuft, die Vorzeichen um bei 


9,9 — 5,4 _—5,3 
+14_+23 +54 
Die Striche bezeichnen diejenigen Bereiche, in denen 
reelle Wurzeln der Gl.(7c) ausgeschlossen sind; sie 
wurden durch Betrachtung der Vorzeichen ermittelt. 
Durch Vergleich der Grenzwertfolgen (12) und (13) 
beschränken sich die Bereiche, in denen reelle Wurzeln 
der Gl. (7) ausgeschlossen sind, auf 
81T 54 5,30 
1AR22735 4,4874, 35 
Die offenen 5 Bereiche sind zwar noch kein Beweis für 
das Vorhandensein von mindestens 5 reellen Wurzeln, 
aber sie machen es wahrscheinlich und erleichtern die 
Suche, die jetzt durch probeweises Einsetzen, z. B. 
nach“dem Horner- Schema, aufgenommen werden 
kann. Man findet dann die 6 reellen Wurzeln 


. (13). 


—6, —4 —2 +1: +3 +5. 
Beispiel2: 
z°6+29°+109°+102+8=0 . (14a). 
Die Umformung gemäß (1b) ergibt 
2 + D=—-10(0 + 5 + 0,8). . . (14b). 


(@+us+p) (+0) = 


wo. a 


a? + = 208 En hat Keine Tesllon Tas gen, ei 


daher keine Rue, u BES) Don nur aus 
. (18) 


j a8 —=0 u MA NE 
und aus TR 
en ı) RAR . (16). 


an 


Für alle x < —1 ist Gl. (14b) schon. voreichenmäße 


unerfüllbar, ebenso für alle > 0. Für —1<x <o 


‚ist Gl. (14b) aber zahlenmäßig nicht erfüllbar. Somit 
ist der strenge Beweis erbracht, daß 
G1. (14) keine reellen Maele hat, Sie > 
‘ wurde aufgebaut aus \ FR 


(0 + 204 ER 


In diesen Beispielen schälten sich die Bereiche aus 


den Umformungen heraus, die nach Belieben vorge- 

nommen werden dürfen. Bei den Gleichungen 3. biesga 
5. Grades kann man sogar einen Bereich vorschreiben. 
Er sei durch r und s begrenzt. Wirbilenr +3= — u 


undrs=p und führen v=a— wein. 
Für die Gleichung 3. Grades gilt dann 
(pPp-+ uwv—b)& + PV—c. 
Für die Gleichung 4. Grades bieten sich 
(++ PR +oVa+g 
=(pvtgu—e)e+pga—d, 


x 


oder =@+g+uv—b) + pad, 
worin g= p+(e—ap)iw 
oder . =alp+gHuu—) at Pu tgu—e), 
worin q=d/p. 
Für die Gleichung 5. Grades gilt 
(a +Uuc + pP) +VC+g) 
= pvrgqu—)R+Pg—dNa—e, 
woring=b—p—uv.. 
Schreibt man hingegen keinen Bereich vor, sondern 


läßt die Bereiche aus der Umformung hervorgehen, 
so hat man bei den Gleichungen 3. bis 5. Grades viele 


Möglichkeiten, von denen nur einige noch Auletuhrt 
seien. 


Umformungen der Gleichung 3. Grades 
(«+ a) (@® +b)=ab—c 


(- aa = —b3—c 
2 +c=—a(ar-+b) 
(a +b)=—an®?—c 

der Gleichung 4. Grades 
A+b®+d=—zlan-+te) 
“+d=— (a0 +bx-+ ec) 
(cc +ta)= —br—ca—d 
(0? + b) (a? + db) = — (au —dx/b+ c) 
(®+b)(® +a)n= (ab—c)u—d, 


der Gleichung 5. Grades / \ 
+) +ar+b)=(ac—d)z+be—e 
(2 +b)(c+a)=(ab—c)2— ds —e 
se Hrbr Hd= —art—ca!—e 
(ra) +bR+d,=(ab—c)® tad—e 
(sta) (#"+ba°+ela) =[ab—c)x + ela—d]x. 


Die sämtlichen Umformungen einer Gleichung lassen 
sich nebeneinander anwenden und ihre ‚‚Grenzwerte‘‘ 
zu einer einzigen Reihe zusammensetzen, durch die die 
reellen Wurzeln oft so stark eingeengt werden, daß: 
man geradezu Näherungswerte erhält, jedoch müssen 
sie jeweils durch Einsetzen geprüft werden. Das be- 
schriebene Verfahren gibt für das Nichtvorhandensein 
reeller Wurzeln in gewissen Bereichen strenge Beweise, 
für das Vorhandensein in den dazwischenliegenden Be- 
reichen zwar nur Wahrscheinlichkeiten, die aber für das 
Aufsuchen der Wurzeln nach weiteren, bekannten Ver- 
fahren wertvollste Dienste leisten und dem praktischen 
Rechner daher willkommen sein dürften. 


Potsdam-Babelsberg. Ei Baula. 
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© NACHRICHTEN 


' chen 1945, 0. Prof. 


p Hochschulnachrichten. : 


In den Jahren- 1945/46 traten an den en 
Universitäten und Hochschulen unter den Dozenten ' 
für Mathematik und theoretische Physik unter anderen 


folgende Veränderungen ein; 


Aachen: Der apl. Prof. a. d. Techn. Hochschule, 


Prof. Dr. R. Sauer zum o. Prof. ernannt. 

Der o. Prof.a.d. Techn. Hochschule Breslau Hubert 
Cremer erhielt einen Lehrauftrag. 

Dr. Josef Meixner wurde zum a. o. Prof. der 
Theoretischen Physik berufen. 

Berlin, Universität: Zu o. Prof. enden aa 


Dr. Hermann Ludwig Schmidt (Mathematik), Dr. 


Kurt Schroeder (Angew. Mathematik), Dr. Robert 
Rompe und Dr. Kurt Möglich (Theoretische Physik). 

Einen Lehrauftrag erhielten: Prof. Dr. E. Rembs 
(Darstellende Geometrie), Studienrat Dr. Schürer 
(Geschichte der Mathematik und Naturwissenschaften). 
- Verstorben ist Prof. Dr. Harald Geppert. 7 

Berlin, Technische Universität: Prof. Dr. 
E. Mohr und Dr. E. Rembs wurden zu o. Professoren 
ernannt. 

Gefallen ist Prof. Dr. W. Cauer. 

Bonn: Im März 1946 verstarb Prof. Dr. 
Bessel-Hagen. 

Braunschweig: Einen Lohraufbrag erhielten: a 
habil. Hilmar Wendt (Mathematik), Prof. Dr. C. 
Weber (Dresden) (Bendetgebie der technischen 
Mechanik). 

Darmstadt: Im September 1944 fielen die Assi- 
stenten Dr. Heinz Thomas und Dr. Ulrich SEROERN" 
witz: 

Dresden: Als Nachfolger von Prof. Dr. C. Mieber 
wurde Dr. H. Neuber zum o. Prof. der Mechanik, als 


Erich 


Nachfolger von Prof. Dr. Ludwig, Prof. Dr. Ött- 


Heinrich Keller zum o, Prof. der Geometrie ernannt. 
Es starben Prof. Dr. M. Lagally im. Januar 1945, 
Prof. Dr. B. Schilling im Sommer 1945, Prof. Dr. 
E. Naetsch im Frühjahr 1946, Prof. Dr. W. Ludwig 
im Dezember 1946. 
Frankfurt: Dr. W. Franz wurde mit der Ver- 
tretung eines Ordinarius beauftragt. 
Fräulein Dr. Ruth Moufang erhielt eine Dozentur. 


Prof. Dr. W. Lorey wurde zum Honorarprofessor in ° 


der sozialwissenschaftlichen Fakultät ernannt. 
Freiburg i. Br.: Dem Dozenten a. d. Universität 
Göttingen Dr. H. Görtler wurde im Dezember 1944 


‚die a. o. Professur für angew. Mathematik übertragen. 


Der a. o. Prof. a. d. Universität Greifswald G. Bol 
erhielt eine Diätendözentur. 
„Prof. Dr. E. Sperner (Straßburg) wurde mit der 
Vertretung eines Ordinarius beauftragt. 

Prof. Dr. Jos. E. Hofmann erhielt: eine Gast- 
professur. 


Gießen: Es starben Honorarprof. Dr. Ph. Maenn- 


‘ Berlin Dr. Günther Schulz wurde als Nachfol er von 
2 - Braun, Dr. H. Wittich, Dr.H. Teer 


 Tollmien (Dresden) berufen. 


 fessor der Mathematik und ihrer Anwendungen RN 


1946,00 0 
Göttingen: Z06 Profe 
Dr. F, Rellich (Dsenden) 
(Königsberg). jr 
Im W.-8. "1946/47 hielt Prof. Dr 
als. Gastprofessor. 7 et 
Zu apl. Professoren. Sie an 


Zum 1. April 1947 wurde Prof. Dr. Ludwig 
emeritiert. Als sein Nachfolger wurde. Prof. ; 


Hamburg: Die Leitung des neugegrü 
stituts für angewandte Methematik wurde ( 
‘H. Zassenhausen übertragen. ; 

Im Spätwinter 1946/47 starb Prof. Dr. =. ne 

Hannover: Rektor ist 2. 2. Prof. Dr. POORER 
Müller. j gut 

Heidelberg: Aus Frnkit wurden die 0. EUR, 4 
soren Dr. H. Seiffert und Dr. W. Threlfall berufen. 

Karlsruhe: Rektor ist z. Z. Prof. Dr. Th. Pöschl. Fe 

Prof. Dr.K. Klotter (T.H. Berlin)wurde zum o. Pro- 


Gefallen ist Dr. Roßbach. .s 
Köln: Umhabiliert hat sich Doz. Dr. M. Pinl. Reha. BT 
bilitiert wurde Prof. Dr. E. Fischer. | 
Mainz: Als o. Professoren wurden berufen: Prof. 
Dr. R. Furch (Rostock) und Prof. Dr. G. Koethe 
(Gießen) als a.o. Professor Dr.H. Wielandt (Tübingen). Re 
Prof. Dr. H. Rohrbach (Prag) wurde eine Gast- N ü 
professur übertragen. ET: 
Prof. Dr.K. Bechert (Gießen) wurde also. Professor RN 
für theoretische Physik berufen. RE 
Es habilitierte sich Dr. F. Oberhettinger: TEE 
Einen Lehrauftrag erhielt Prof. Dr. C. Schmieden KERN‘. 
(Darmstadt). er 


Marburg: Dr. H. Grötzsch habilitierte sich von x 
Neuem. Prof. Dr. E.R. Neumann wurde emeritiert. Ert: 


München, Techn. Hochschule. Einen Lehrauf- 
trag erhielt Dr. W. Damköhler. 

Münsteri. W.: Berufen wurde der o. Professor a. d. 
Universität Jena Dr. F.K. Schmidt. Er wird als 
Gastprofessor Vorlesungen a. d. Universität Berlin 
halten. Mit seiner Vertretung wurde der Doz. Dr. x 
Th. Schneider (Göttingen) beauftragt. 

.Dr.'H. Ulm wurde zum a. o. Professor ernannt. 

Prof. Dr. H. Behnke ist für ein Semester als Gast > 
an die Universität Upsala eingeladen. Rt: 

Rostock: Der o. Prof. der angew. Mathematik an 
der Universität Halle Dr. Harry Schmidt erhielt 
einen Ruf auf das Ordinariat der Mathematik. 

Stuttgart: Als Nachfolger von Prof. Dr. W. Kutta 
wurde Prof. Dr. F. Lösch (Rostock) berufen. 

Dr. O0. Baier wurde zum a. o. Professor ernannt. 

Tübingen: Als o. Professor wurde Prof. Dr. E. 

Kamke rehabilitiert. Prof. Dr. G. Hamel (Berlin) 
wurde eine Gastprofessur übertragen. > 


Betr.: ERNEUERUNG DES ABONNEMENTS 


Wir bitten unsere Bezieher an die rechtzeitige Erneuerung des Abonnements zu denken, 
soweit zunächst nur ein Vierteljahres-Bezug bestellt wurde. \ 


} 


Bei Post-Abonnement muß die Erneuerung bis spätestens zum 26. des letzten Quartals- 
Monats bei dem zuständigen Postamt erfolgen. Bei Bezug durch eine Buchhandlung oder 
durch den Verlag bitten wir den Buchhändler oder den Verlag rechtzeitig zu verständigen, 
daß das Abonnement verlängert werden soll. 


nn 


Verantwortlich für den Inhalt: Prof. Dr. Fr. A. Willers, 


Dresden; ‘für den Verlag: H. Kaesser, Berlin. 


Verlag: Akademie-Verlag GmbH., Berlin N 4, Chausseestr. 106, Fernsprecher: 42 50 01 (Verlag App. 274, Vertrieb 


App.275), Postscheckkonto: Berlin 35 021. 


Bestell- und Verlagsnummer dieses Heftes: 1009/25— 27/3. 


Die Zeitschrift 


für angewandte Mathematik und Mechanik erscheint monatlich. Bezugspreis: Vierteljährlich RM 15.— zuzüglich Bestell- 


geld, Einzelheft RM 6.—. 


Verantwortlich für den Anzeigenteil: 


Curt F. W. Schreiber, Berlin. Druck: 


Julius Beltz, Langensalza. Veröffentlichtunter derLizenz-Nr.245 der Sowjetischen Militär- "I OTwWaltung in Deutschland. 


(Th. B.-Nr. 04011/16 b) Julius Beltz, Buchauue 2, Langensalza. — Nr. 888. 


L.-Nr. 245 


